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XII. The Waveforms of Atmospherics and the Propagation of very 
Low Frequency Radio Waves. 


By P. W. A. Bowe, 
Cavendish Laboratory, Cambridge *. 


[Received November 28, 1950.] 
[Plate VI.] 


SUMMARY. 


The responses of narrow-band receivers to individual radio atmospherics 
have been observed in order to describe the propagation of radio waves 
of frequency 2 ke./s. to 10 ke./s. over the surface of the earth. It has been 
found that the atmospherics from a fixed distance are sufficiently 
uniform to enable the relative attenuations of different frequencies 
to be deduced. The study reveals that waves of frequency below about 
8 ke./s. are heavily attenuated during the daytime but are propagated 
freely at night. The attenuation is greater during a sudden ionospheric 
‘disturbance. 

These results confirm and extend those obtained recently by Gardner. 
‘They can be explained in outline on a theory recently proposed by 
Budden. 


$1. INTRODUCTION. 


SEVERAL workers have shown that the waveforms of some atmospherics 
can be interpreted in terms of irregularly shaped radio impulses reflected 
one, two and more times from the ionosphere (Laby, McNeill, Nicholls, 
Nickson, 1940 ; Schonland, Elder, Hodges, Philips and van Wyk, 1940). 
Under good conditions it is often possible to deduce the height of 
reflection of these impulses, and heights of about 85 km. have been 
deduced. 

_Other workers (Bracewell, Budden, Ratcliffe, Straker and Weekes 
1951) have investigated the propagation of very long radio waves emitted 
by commercial senders, and a considerable body of knowledge now 
exists concerning the height from which these waves are reflected by 
the ionosphere, and their intensity. 

It is clearly of interest to compare the results of the two types of 
experiment. Unfortunately, when this comparison is attempted, a 
serious difficulty arises because in the atmospherics observed the shape 
of the impulse which is successively reflected from the ionosphere is 
usually irregular, and varies from sample to sample, so that it is impossible 
to specify with any accuracy the radio frequency to which it corresponds. 


Se 
* Communicated by Mr. J. A. Ratcliffe. 
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In the work to be described in this paper radio atmospherics are used 
in a different way to investigate the propagation of very long radio. 
waves. A single atmospheric is applied simultaneously to four narrow- 
band receivers tuned to different very low frequencies between 2 ke./s. and 
10 ke./s., and the responses are compared. The relative responses on the 
different frequencies (which might be called the spectrum at the receiver). 
will depend on (a) the relative intensities excited in the different 
frequencies by the particular lightning flash being measured (which. 
might be called the spectrum at the source), and (b) the characteristics 
of the transmission path between the flash and the receiver. In order 
to separate these two factors, the atmospherics received from flashes 
at different distances have been compared. 

The assumption is made that, statistically, all flashes excite the same 
spectrum at the source, at least in the region of very low frequencies 
here considered, and a search is made for some significant change 
produced in the received spectrum as the distance of the flash changes.. 
For one of the two classes of atmospherics received, described in more 
detail below, it has been found possible to establish a change of this kind 
accompanying a change of distance. This provides confirmation of the 
original assumption that there is some statistical constancy in the 
spectrum produced at the source by different lightning flashes. That 
change in the spectrum which accompanies increase of distance is 
attributed to the propagation conditions. It is also found that the 
propagation conditions, deduced in this way, change from day to night 
and during a sudden ionospheric disturbance. The results are summarized. 
and discussed in § 5. 

The recorded atmospherics fell into two easily distinguishable classes 
as follows. 

Class I, had a fairly smooth “waveform” with marked low frequency 
periodicity of a more or less regular nature. The narrow-band amplifiers. 
tuned to frequencies below 10 ke./s. gave a large response with. 
atmospherics of this class. An example is shown in Pl. VI. (i.). 

Class IT. had a highly irregular “‘“waveform’’ which showed an erratic: 
variation with many short period fluctuations. When the overall 
amplitude of this type of atmospheric, as recorded on the wide-band 
receiver, was of the same order as that of an atmospheric of Class I., the 
response of the narrow-band amplifiers was much smaller for frequencies 


below 10 ke./s. An example of Class II. is shown in PI. VI. (ii.) and can. 


be compared with the atmospheric of Class I. shown in Pl. VI. (i.) 


It appears that these two classes of atmospheric are different at the 
source, and in this paper only those of Class I., which gave an appreciable- 
response on the narrow-band receivers, are considered. They formed. 


about 75 per cent of the atmospherics recorded. The two different 


classes have been noticed by most workers in the subject, and Appleton. 


and Chapman (1937) have ascribed them to the leader (Class IT.) and 


return (Class I.) strokes of the lightning flash. For the purposes of the: 


t 
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present paper, it is sufficient to distinguish between them by the nature 
of the recorded waveforms. 

It will be shown in §4 that, when atmospherics of Class I. are selected 
in this way, their spectrum at the source is, within limits, the same for 
all atmospherics, so that changes in the spectrum due to propagation. 
conditions can be shown up. Since atmospherics of Class IT. produce: 
only small radiation at the frequencies here considered, we assume that 
they do not radiate appreciably at these frequencies and therefore have 
a different spectrum at the source. Both classes of atmospheric have been 
received from distances up to 1700 kilometres.’ 

The sources of atmospherics which were received and analysed in 
the way here described were simultaneously located by radio direction 
finding methods by the “ Sferic” section of the Meteorological Office. 
This provided the measure of distance to the source of the atmospheric. 

In the work described here each atmospheric was received on four 
narrow-band receivers tuned to four different very low frequencies, and 
for the sake of comparison, the waveform of the atmospheric as received. 
on an aperiodic amplifier was simultaneously recorded. The responses 
at the four spot frequencies represent four points on the spectrum of 
the received atmospheric and could be deduced if that spectrum were 
known. In several cases the spectrum of the received atmospheric has. 
therefore been deduced from the waveform by numerical Fourier 
analysis and it has been shown that the responses on the four tuned. 
receivers correctly represented the four points on the spectrum. It 
was at one time thought that the method of numerical Fourier analysis. 
of the waveform would be simpler than the method of reception on four 
narrow-band receivers, but it has proved that the latter method is in 
fact simpler and quicker. Since a large number of results are required 
for making a statistical analysis, the method of observation with the 
four receivers has therefore been used. 

The results described in this paper confirm and extend those recently 
obtained by Gardner (1950). He recorded the integrated “level” of 
atmospherics received on four narrow-band receivers tuned to four 
different very low frequencies and observed how these levels altered with 
time of day and with the distance of a storm centre. He deduced 
information about the propagation conditions on the different frequencies. 

While he was working in Cambridge, Garduer also constructed the 
equipment, described in outline in this paper, for recording the response 
of narrow-band receivers to individual atmospherics, and obtained some 
preliminary results with it before the present author took over the work 
from him. It is hoped that the equipment will be described in detail 
i Gardner. 
si apis given in this paper, and the corresponding results of 
Gardner, have been explained, at least in outline, ihe Budden (1951). 
A preliminary comparison of his theory with the present results is made: 
in §5. ; va tage 
re K 2 
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§ 2. EXPERIMENTAL METHOD. 


The atmospherics were received on a vertical aerial and passed. through 
a cathode follower to an aperiodic receiver and four tuned receivers all 
in parallel, as shown in fig. 1. The output of each receiver was applied 
to a separate cathode ray tube and when an atmospheric arrived all 
five tubes were photographed. Thus the waveform of the transient 
was recorded together with the response of four tuned receivers to it. 


Fig. 1. 


BRIGHTENING 


CAT HODE 
FOLLOWER 


AMPLIFIER 


To other Tuned 
Amplifiers with 
C.R.0. Displays 


HIGH VOLTAGE SAWTOOTH 
SUPPLY TIMEBASE 


Diagram showing the experimental arrangement for recording atmospherics 
on an aperiodic and several narrow-band receivers. 


2.1. Receiving Unit. 


The vertical aerial, which was erected on a site near Cambridge 
relatively free from noise, had a length of 9 metres and a capacity to 
ground of 120 puF. The voltage in the aerial was applied to the grid 
of a cathode follower across a grid resistor of 10° ohms, situated at the 
bottom of the aerial. A long cable connected the low impedance output 
of the cathode follower to a pre-amplifier with a step attenuator, the 
output of which was applied to the five receivers in parallel. 

The aperiodic receiver had a gain which was constant to within 1 db. 
for frequencies between 800 c./s. and 50 ke./s, and gave an output which 
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was proportional to the electric field at the aerial. It had two stages of 
amplification followed by cathode and anode followers to give a push-pull 
output to the cathode ray tube. 

Each tuned receiver had a bandwidth of about 180 c./s..between half- 
power points, and each was tuned to a different frequency somewhere 
between 2 ke./s. and 10 ke./s. Each receiver consisted of an amplifier 
with a tuned-grid stage followed by a second amplifier, from which 
positive and negative feedback were applied to the first amplifier in 
order to narrow the bandwidth. Each receiver was followed by anode 
and cathode followers to give push-pull output for the cathode ray tube. 


2.2. Display Unit. 


The output from each of the four receivers was applied to the 
Y deflection plates of four 3-in. cathode ray tubes. The same sawtooth 
time-base voltage was applied to the X plates of each tube ; the time of 
sweep was about 2 milliseconds and the time of fly-back was less than 
5 microseconds. This time-base was freely running so that the atmospheric 
appeared at a random place on the trace. 

The grids of the cathode ray tubes were normally maintained at 
a negative potential so that no trace appeared. The output of the 
aperiodic receiver was connected to a brightening unit which, when an 
atmospheric above a certain small amplitude arrived, applied a positive 
pulse to the grids so as to brighten the tubes for the duration of one time- 
base sweep, and this trace was automatically photographed. The film 
was moved on manually and, to avoid double exposure, a delay was 
provided in the brightening unit so that, after brightening once, the 
circuit was inoperative for about half a second. 

A typical record obtained in this way is shown in Pt. VI. (i.). The 
waveform as received by the aperiodic amplifier is shown on the centre 
trace and the responses of the four narrow-band receivers on the other 
traces. Each trace brightened at a point (a) and continued bright to 
(b), whence the fly-back returned the spot to (c) and the rest of the 
time-base was completed. 


2.3. Measurement of Distance. 


The distance to the source of the atmospherics was determined by 
taking advantage of the “ Sferic ” direction finding organization of the 
Air Ministry Meteorological Office (Ockenden 1947), and the grateful 
thanks of the author are due to the Director and staff for helpful 
co-operation. 

Use was made of the routine service by which the sources of individual 
atmospherics are located by the intersections of radio bearings taken 
simultaneously at four widely-spaced observing stations. The use of 
telephone connection during an experiment enabled individual 
atmospherics recorded at Cambridge to be identified with sources located 
by the “ Sferic ” organization. 
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The distance to the source was then measured directly. There was 
‘a small chance of mis-identification, but this was not serious enough to 
affect the results. The sources could generally be located to an accuracy 
ot about 10 per cent up to 1000 km. and less accurately beyond, but 
this was quite sufficient for the present purpose. 


§ 3. Tae SPECTRUM OF THE RECEIVED ATMOSPHERICS. 


We shall here discuss the spectrum of the received atmospheric as 
deduced from Fourier analysis of the recorded waveform and relate it 
to the observed responses from the narrow-band receivers. 


3.1. The Output from the Tuned Receivers. 


It is interesting to notice the way in which the narrow-band receivers 
respond to the input transient depicted at (2), Pl. VI. (i.). If we look, for 
example, at the recorded response of the receiver tuned to 7-5 ke./s. 
we see that, while the original transient is arriving, the amplitude is. 
building up so that, when the transient is over, the circuit is left 


performing a decaying oscillation with a decay time determined by the ~ 


band-width of the amplifier. The amplitude of this oscillation represents 
the amplitude of the appropriate Fourier component in the analysis 
of the transient. The amplitude of the final decaying oscillation cannot, 
-of course, be decided until the transient is over, so that the response 
‘during the building-up period need not bear any relation to the final 
amplitude when the transient is over. The record of the response on 
3:5 ke./s. shows this well: during the “‘ build-up” period there is quite 
.a large response, but when the transient is over, the circuit has settled 
‘down to give a negligibly small decaying oscillation, corresponding to 
the fact that the Fourier component of frequency 3-5 ke./s. is too small 
to record. 

The four narrow-band receivers were calibrated by observing their 
responses to a test signal which consisted of a rectangular pulse of voltage 
of known amplitude and of duration 14 psec. This pulse was repeated 
fifty times per second but, of course, only occasional ones were recorded. 
Each of these pulses had an effectively uniform spectrum over the 
relevant range of frequencies 2 ke./s. to 10 ke./s. Although it is considered 
that the absolute response of the different receivers could be made equal 
by this method of calibration, it will become apparent, in what follows, 
that the conclusions to be reached in this paper do not depend on this 
adjustment for equality. All the measurements were made in terms of 
the ratio of the responses at two frequencies, and the way in which that 
ratio changes from time to time. It is only necessary, then, that the 
ratio of the responses should be determined, and this was done by the 
method of calibration described. 

A similar argument helps us to avoid a difficulty in the measurement 
of the responses on the narrow-band receivers. There might be some 
doubt as to whether the amplitude of the decaying oscillation should be 
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measured at the time when the transient is finished, or whether the decay 
should be allowed for, and the amplitude be estimated at the time when 
the transient started. Since, however, we are only concerned with the 
ratio of amplitudes on the different frequencies, it does not matter which 
Measurement is taken, provided it is kept the same for all frequencies. 

Records of the type illustrated in Pl. VI. were made on photographic 
film and were measured by projecting and comparing with a scale. 
3.2. Direct Fourier Analysis of the Transient. 

The method described in the previous section could be thought of as 
giving four points on the frequency spectrum of the received atmospheric. 
It was at one time thought that this information could be more easily 
obtained by performing a Fourier analysis of the transient recorded on 
the wide-band amplifier, and, moreover, it appeared that if the analysis 
were performed in this way, attention need not be restricted to four 
frequencies. Some of the records of transient responses were therefore 
analysed by ordinary methods of Fourier analysis, as described below. 
It soon became apparent, however, that the labour of performing a 
full Fourier analysis was considerable, and, since the aim of the work 
was to record numerous atmospherics from which statistical deductions 
could be made, this method was not continued. It proved to be much 
simpler and quicker to use the apparently more complicated method 
of the four tuned receivers to obtain the Fourier components at the four 
chosen frequencies. 

It is, however, of interest to see how far the Fourier analysis of the 
recorded wave-form agreed with the magnitudes of the Fourier com- 
ponents deduced from the responses of the four narrow-band amplifiers. 
For this purpose the transient waveform was redrawn on as large a scale 
as possible and three standard methods of Fourier analysis were used. 
‘These were :— 

(1) Anelectro-mechanical analyser developed by the Admiralty Research 

Laboratory (Barber, Ursell, Darbyshire and Tucker 1946). 
(2) Beevers Lipson Strips (Ross 1943). 
(3) Hollerith Computing Machine (Hodgson, Clews and Cochran 
. 1949). 

Fig. 2 (a) gives the results of applying methods (1) and (2) to the 
waveform shown and fig. 2 (b) shows the results of applying methods 
(2) and (3) to another waveform. It will be seen that in each case the 
deductions made by the two different methods agree quite well. The 
points marked by crosses show the relative magnitudes of the Fourier 
components as determined by measurement of the responses on the 
narrow band receiver. These magnitudes have been adjusted so that 
the point for 7:5 ke./s. falls on the computed curve. It will be seen that 
they agree well with the results calculated by direct Fourier analysis 
of the transients, and we are satisfied that the method of using the four 
tuned receivers gives reliable results for the ratio of Fourier components 


at the relevant frequencies. 
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Fig. 2 (a). 
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received spectra of two atmospherics as deduced by the methods indicated.. 


The four points deduced from the responses of the tuned receivers are 


indicated by crosses. The waveform analysed is also shown. 
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It might be thought that an important part of the atmospheric wave- 
form might fail to be recorded on the wide-band receiver because of 
the time lost before the cathode-ray trace was brightened by the 
atmospheric. The good agreement between the spectra computed from 
the recorded waveform and the responses of the four narrow-band 


receivers shows that there is no difficulty of this kind in the range of 
frequencies here considered. 


§ 4. ANALysts or RESULTS. 
4.1. Method of Analysis. 


In order to clarify the method of analysis, it is convenient to define. 
some attenuation factors. Consider a lightning flash which has a radiation 
field Ad(f) at the source, where A depends on the strength ot the source. 
and ¢(f) represents the way in which the field varies with the frequency /f. 
Then in free space the field received at a distance r would be given by 
A¢(f)/r. Suppose now that this atmospheric were transmitted a distance 
r over the surface of the earth and that as a result each component 
frequency f suffered an attenuation p(f, 7); so that the received field would 
be given by 


Bf. == 8) pln 


Now let the field produced by a source at a distance r be measured. 
at a series of different frequencies f,, /;, fo . . . f,, and let the result be 


expressed as a ratio between the field at the frequency f,, and the field. 
at a reference frequency f, so that we obtain the quantity 


Ef, r)/E(fo; 1) ={b(fnr)/ (fo) $0 Tne r)/p(fos ry}. 


Next, suppose that these measurements are made for several atmospherics. 
at a series of different distances, and that the ratio E(f,,, 7)/E(fo, 7) is 
plotted as a function of r. The resulting curve will show how the ratio. 
p( fr» *)/p(fo. 7) of the attenuation factors at the two frequencies Sn and fp 
varies with distance. Since the attenuation factor p(/, 7) is unity for all 
frequencies when the distance is small, the intercept on the axis at 
r=0 will give the magnitude of 4(f,,)/¢(/o), Which depends on the spectrum 
of the atmospheric at the source. 

It is of course most unlikely that the spectrum 4(f) of the atmospheric- 
at the source is the same for all atmospherics. If, however, they do not 
differ too much among themselves then, by taking enough observations | 
on atmospherics originating at all distances, it may be possible to draw 
a smooth curve through the set of observed points so as to show, with 
some real significance, how the ratio p(f,,7)/p( To: r) of the attenuation 
factors varies with distance. We shall show that, in fact, the BpaeIT 
d(f) are sufficiently nearly the same for frequencies between 2 and 
10 ke./s. for this procedure to be of use. 
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4.2. Measurements under Normal Conditions. 

In order to examine seasonal changes the observations are grouped 
in four sets as summarized in Table I. which lists the frequencies used, 
the number of experimental runs, the number of atmospherics recorded, 
and the distances between which the sources were located. 


TABLE I, 


- Frequencies No. of No. of Distance of 
Occasion ko./s: 


runs | atmospherics | sources, km. 


Winter day (3-5), 5, 7-5, 10 11 600-3600 
Night 3°5, 5, 7:5, 10 10 BE 1000-3000 
Summer day | 3-5, 5, 7:5, 10 300-2000 


| Summer day | 2:25, 4:25,6-25,] 256 0-3000 
7-5 


In order to analyse the results the frequency 7-5 ke./s. was taken as the 
standard frequency fp, and for each atmospheric observed the ratio 
K(f, 7)/E(fo, 7) was deduced and plotted against the measured distance 


Fig. 3 (a). 


0 400 800 1200 1600 2000 2400 
Distance, km. 


To show how the ratio E(f,, 7)/E(fo, 7) of field strengths at 4-25 ke,/s. (=f,) and 
7-5 ke./s, (=f) for atmospherics at different distances, varies with distance. 
All observations were made in daytime. 


r of the source. A series of plots of this kind was made for each of the 
frequencies observed. Fig. 3 (a) shows the results obtained on summer 
-days for f, = 4-25 ke./s. and fy =7°5 ke./s. 
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It will be seen that although there is a considerable spread in the 
observed points, of the order of about 10:1 at all distances, it is 
nevertheless possible to draw a smooth curve through them which 
clearly has some significance and which shows that the ratio decreases 
as the distance increases. We deduce from this curve that the intensity 
of waves of frequency 4:25 ke./s. decreases more rapidly with distance 
than that of waves of frequency 7-5 ke./s. 

There is, of course, a very considerable spread amongst the observed 
points plotted in the figure. This is caused largely by the differences 
in the ratio 4(f,)/4(fo) from one atmospheric to another, but some of it 


Fig. 3 (0). 
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+ Observations at night. 
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‘is also probably due to inaccuracies in the measurements of the a 
of the source, and to wrong identifications between some o : : 
-atmospherics whose distance was measured and those Pie spec ee 
‘response was determined. Examination of the results on a lees 
and at all times has led us to ae ae we poe a a fo) varies by 
r of about 10:1 from atmospheric 0 atmosp : 

ae 3 (b) in which the ratio E(fy, r)/E(fo/) is ee ee = 
fo=3'5 ke./s. and fy=7-5 ke./s., shows similar Sain \apeae 
iced by day are shown as dots, and a smooth ous - ia S we 
to lie evenly amongst them. Observations taken by night are sho 
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by crosses. It will be seen that, during the day, for distances greater 
than about 1200 km. the response on 3-5 ke./s. is in nearly all cases so 
small as to be negligible and is less than 1/100th of the response on a 
frequency of 7:5 ke.js.* It is clear that, whereas waves of frequency 
7:5 ke./s. can be transmitted over these distances during the day with 
comparative ease, those of frequency 3-5 ke./s. are heavily attenuated. 
This result, shown clearly by these observations alone, is the most 
striking feature of the transmission characteristic which we shall 
summarize in the next section. 

The other important point which arises from fig. 3(b) is concerned with 
the marked difference between the results for day and night. It is clear 
that although waves of frequency 3-5 ke./s. cannot be transmitted over 
distances greater than 1200 km. by day they are transmitted almost as 
easily as those on 7:5 ke./s. by night. — 


Fig. 4. 


E(f.7) 
EG.) 


F=7 5k 
I 


0 400 800 1200 1600 2000 2400 
Distance, km. 


Curves showing how the average ratio of field strength E(f,,. r)/E(f,, 7) varies 
with distance for fy=7-5 ke./s. and f, as ee ete 


Results obtained with the other frequencies have been plotted in the 
same way and smooth curves have been drawn through the points. 
The complete plots are not reproduced, but the smooth curves are shown 
in fig. 4 for the summer day observations. These observations were 
nearly all made within two hours of local noon, and most of the 
atmospherics came from directions between 030° and 135° EB. of N. The 
analysis was based on those recorded atmospherics for which a definite 
location was obtained by the Sferics organization. 

The night and winter-day observations did not include a sufficient 
number of atmospherics to enable similar curves to be drawn, and they 


ES EASE yon et a Polk Su eT i BOE 

ae The rena Bare on the two frequencies of 7:5 and 3-5 ke tone 
istant atmospheric in the daytime is well illustrated by th jen 

In PRSY tli y. y the example shown 
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are not presented here in detail. However, the observations on 
3-5 ke./s. already mentioned do establish the important fact that the 
received amplitude E(f,,7r)/E(f,7r) for frx=3'5 ke./s., fo=7°5 ke./s., 
r=1000 km. is at night about ten times what it is during the day. This 
means that the attenuation factor p(f2, r)/p( fo, 7) is at night ten times what 
it is during the day. has 
4.3. Deduction of Attenuation Curves. 
The curves of figs. 3 and 4 show how the quantities 


(P(Sn)/b(fo)} Lota» ")/p( fo, ”)} 


vary with distance when f, is taken as 7-5 ke./s. It is now desirable 
to normalize these curves so that they pass through the value unity for 
small distances; they then give information about the attenuation 


Fig. 5. 
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Curves showing how the ratio of attenuation coefficients p(f,,7)/p(fo, 7) varies 
with distance for fy)=7:5 kc./s. and f, as indicated. 


coefficients alone, but not about the ratios 4(f,,)/d(fo). The normalizing 
factor was obtained by extrapolating the curves to uae pian Sa 
It is also convenient to take a new frequency, 10 ke./s., for /o*, and we 
then obtain, for summer daytime, the series of curves shown in fig. 5 in 
which p(f,, 7)/p(fo, 7) is plotted against distance with Jo=10 ke./s. 

It is convenient to present the same results in still another way, He 
shown in fig. 6, in which the ratios p(f, r)/p(fo, 7) are ST Resin ets 
frequency f for three different distances r= 500, 1000 an s a 
In these curves f, was taken as 10 ke./s. and if there were no a nee 
at this frequency, the curves would represent attenuation as a function 


i i rs i , F- 
* In the preliminary calculations it was more convenient to take fy=7'5 oe s 
because one of the receivers was working on this frequency im every experiment, 


but there was not always a receiver tuned to 10 ke./s. 
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of frequency at the three distances ; if there be some attenuation on 
10 ke./s., the curves showing total attenuation as a function of frequency 
would have their scales altered by different amounts so that they would 
no longer all pass through the same point at 10 ke./s. 

By a different experimental method Gardner (1950) has plotted similar 
curves, and the cross at 5 ke./s. (fig. 6) represents his results for 
atmospherics originating from a distance of about 1500 km. It agrees 
well with the present results. 

Neither the winter nor the night observations are sufficient to enable 
similar attenuation curves to be drawn. The night observations on. 
3-5 ke./s. shown in fig. 3(b) do, however, indicate, as has been mentioned 
before, that the attenuation factor p(f2,7)/p(fo.7) for fs=3°5 ke./s., 


Fig. 6. 
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Summary of attenuation coefficients for different conditions. 

Variation of ratio of attenuation coefficients p(f,7)/p(fo, 7) with 
frequency f for fp>=10 ke./s. and 7 as indicated. Normal day 
conditions. 

---- Similar curve for r=2000 km. during an §.I.D. 
+ Value obtained by Gardner for atmospherics from about 1500 km 
©. Night attenuation on 3-5 ke./s. for atmospherics from 1000 km. 


fo=7'5 ke./s., r=1000 km. during the night is about ten times what it is 
during the day. Since the day attenuation is about 0-05, this gives a 
value of 0-5 for the attenuation at 3-5 ke./s. This value is shown by the 
open circle on fig. 6. 


4.4, Measurements during a Sudden Ionospheric Disturbance. 


It has been noted by Gardner (1950) that, during a Sudden Tonospheric 
Disturbance (S.I.D.), the integrated level of atmospherics decreases on 
frequencies below 7 kc./s.,so that the spectrum of anindividual atmospheric 


——— 


a 
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might be expected to show a similar change. Fortunately it was 
possible to make observations during a large S.I.D. on 20th February 
1950. Atmospherics were recorded at five separate times before, during 
and after the 8.I.D., as indicated in fig. 7. The set of observations at 
1550 G.M.T. was made while the “ Sferic ” organization was operating, 


Biga7: 
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a) Change in ratio of received field strength of atmospherics, 


E(fa, 7)/E(fo, 7) when 8.I.D. occurs. 
b) Change in ratio of received field strength of atmospherics. 


E(f,, 7)/E(fo, 7) when 8.I.D. occurs. 
(c) Sudden enhancement of integrated level of atmospherics on. 


22 ke./s. 
(d) Beginning and end of fadeout reported by B.B.C. 


Crossss indicate mean values for sets of points. 
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so that definite distances were available for the atmospherics recorded. 
In addition the routine ‘“Sferic” charts of the Meteorological Office 
indicated that there was a strong source of atmospherics at a distance of 
2400 km. throughout the whole day and very little activity elsewhere. 
From this it can be stated, with some certainty, that the recorded 
atmospherics originated in this distant storm. 

The results are shown in fig. 7. In (6) the quantity E(f,)/E(fo) is 
plotted for f,=5 ke./s. and fy=10 ke./s., in (a) the same ratio is plotted 
for f,;=7:5 ke./s. and fy=10 ke./s., in (c) the integrated “level” of 
atmospherics as recorded on 22 ke./s. is shown, and (d) shows the extent 
.of the “ fade out ”’ as reported by the B.B.C. In (a) and (6) each recorded 
atmospheric is indicated by a dot, and the mean values for the set of 
dots in each series is shown by a cross. The curves are drawn through 
the crosses. The most striking effect is seen in curve (b) which shows 
‘that, near the peak of the S.I.D., the component of frequency 5 ke./s. 
decreased about ten times compared with the component of frequency 
10 ke./s. This is in good agreement with Gardner’s observation that the 
integrated “level” of atmospherics as received at 5 ke./s. may often 
become too low to observe during an §.I.D. Although the corresponding 
‘decrease on 7:5kc./s., shown in curve (a), is also noticeable, it is not so 
marked. 

It is possible to apply the results shown in fig. 7 to the daytime 
-attenuation curve of fig. 6 so as to derive an attenuation curve which 
would represent the conditions during an intense 8.I.D. Since fig. 7 applies 
to a distance of about 2400 km., we use it to modify the attenuation 
curve of fig. 6 drawn for the approximate distance of 2000 km., and we 


-arrive at the dashed curve in fig. 6 to represent the conditions during 
the 8.I.D. 


§ 5. CONCLUSIONS. 


The final conclusions of this work are represented by the attenuation 
‘curves of fig. 6. These curves show how the attenuation factor p(f, 7) at 
different distances varies with frequency relative to the attenuation 
p(fo, 7) at a standard frequency f,. The three full curves are for normal 
‘daytime conditions, the dotted curve for abnormal daytime conditions 
‘during a sudden ionospheric disturbance and the open circle indicates 
the result for 3-5 ke./s. during the night. 

From the nature of the results it cannot be claimed that the curves 
-are accurate to more than about 50 per cent. It can be seen, however, 
that an accuracy of this order is adequate to establish the general nature 
-of the propagation characteristics of the space between the earth and the 
ionosphere in the frequency range 2 ke./s. to 10 ke./s. 

A first attempt at a theoretical explanation of these results has been 
made by Budden (1951). From consideration of the waveguide formed 
by the earth and ionosphere, he has deduced theoretical curves similar 
ito those of fig. 6. Budden proceeds on the assumption that the 
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ionosphere has a sharp lower boundary at high h, with appropriat 

values of electron density and collision frequency. From Retrem : 
of his theory, attenuation curves for 1000 km. path for heights of 85 iat 
and 50 km. have been obtained. These are shown in fig. 8, to ie 
with the experimental curve for the attenuation 1000 an fee ae 
source. It can be seen that there is semi-quantitative essen! ae a 
the theory and the experimental results. According to Budden’s theor 

the curve will be moved to the right for decrease in reflection hei ht 
and to the left for increase. It is probable that movements of this kind 


explain the different experimental curves obtai : 
and during the night. obtained during an §.J.D. 


Fig. 8. 
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CAPTION FOR PLATE VI. 


Two examples of the record made on individual atmospherics. For both 
examples (i.) and (ii.) the traces are the responses of different amplifiers 
as follows :— 

1. Tuned to 10 ke./s. 

2. Wide-band. 

3. Tuned to 3-5 ke./s. 

4. Tuned to 7-5 ke./s. 
5. Tuned to 5 ke./s. 

Tn example (i.), trace 2, (a) indicates the beginning of the atmospheric when the 
trace is brightened, and (b) is when the time-base starts the fast return. 
to (c). The trace is bright for the duration of one timebase stroke so that 
the trace finishes just before (a). All traces in both (i.) and (ii.) are obtained 
in the same way. 

ery ) is typical of Class I. mentioned in the text and example (ii.) of 
lass IT. 


P. W. A. BOWE 
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SUMMARY. 

‘A calculation is made of the elastic constants of aluminium which 
correspond to pure shears without change of volume, and an explanation. 
is given for the fact that aluminium crystals are almost elastically isotropic.. 
Conclusions are drawn about the nature of the density of states near the. 
Fermi surface, and it is shown how the elastic constants may be expected. 
to depend on electron concentration in aluminium solid solutions. 


§1. INTRODUCTION. 
It is well known that the total energy of a metal may be expressed approxi- 
mately as a sum of separate terms. If we can calculate how these terms 
change as the metal is sheared, we can obtain the elastic constants of the 
metal as a corresponding sum of separate contributions. 

Calculations of the shear constants (7. e. the elastic constants which relate 
to homogeneous strains without change of volume) of single crystals of 
metals have hitherto only been carried out for monovalent metals, namely 
the alkali metals and copper (Fuchs 1936). In these metals two important 
contributions appear, which arise respectively from (i.) the purely electro- 
static correction term which has to be added to the potential energy 
obtained by using the sphere approximation of Wigner and Seitz, and (ii.) 
the term in the energy which is due to the non-Coulomb forces between the 
ion-cores. The first of these contributions is very anisotropic, the anisotropy 
ratio {C44/$(C,,—Cyj2) }a1, s, being 8-9 for a face-centred cubic, and 7-4 for a 
body-centred cubic structure. For these monovalent metals the Fermi 
surface is approximately a sphere lying wholly inside the first Brillouin 
zone. The effect on the Fermi energy of movements of the bounding 
planes of the zone as the metal is sheared is therefore small, and the 
corresponding contribution to the elastic constants is negligible. In the 
case of metals of higher valency, however, this is no longer the case. In 
this paper we make an estimate of the contribution to the elastic constants 
from the Fermi energy in the case of aluminium, this metal being chosen. 
for the following reason. Jones (1949) has pointed out that in the case of 
aluminium the contribution to the elastic constants from the non-Coulomb 
interaction between the ions must be negligible because of the large clearance 
between the ion-cores. The experimental values of the constants for 
aluminium, however, which are almost isotropic, are not consistent with 
their arising solely from the very anisotropic electrostatic contribution. 


* Communicated by Professor H. Jones. 
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We are therefore led to look for another important contribution. The 
room temperature values for aluminium, as given by Lazarus (1949), are 
4(C,,—Cy2)=2-085. 1012 dynes cm.~? and C,,=2°853 . 1014 dynes cm.-?, 

It is shown in this paper that it is in fact possible to account for each of 
these observed values as a sum of contributions from the Fermi energy and 
from purely electrostaticenergy terms. The large positive contribution to C4, 
from the electrostatic terms is balanced, apparently accidentally, by a 
correspondingly large negative contribution from the Fermi energy. In 
order to obtain agreement with the observed values of both C,, and 
4(C,,—C,.), it is necessary to assume that there are some overlapping 
electrons over the square faces of the Brillouin zone as well as over the 
hexagonal faces. ) 

Following on these calculations it is also possible to see how the elastic 
constants would be expected to change with electron concentration as 
metals of higher or lower valency are dissolved in aluminium. There does 
not appear to be any experimental data relating to this. 


§ 2. Tor ToTaL ENERGY IN THE HARTREE APPROXIMATION. 


In metals with more than one valence electron per atom, the approximate 
form for the wave functions of the higher states used by Wigner and Seitz 
for monovalent metals, viz., %o(r) exp(tk.r), where %)(r) is the wave 
function of the lowest state, can no longer be used. It follows that the 
Fermi energy now contains a potential as well as a kinetic part, and the 
valence electron density is no longer given by —e | |?.__ In view of this, it 
is necessary to reconsider the theory of the division of the total energy into 
separate terms. 

In what follows we shall, as a first approximation, neglect valence 
electron exchange and correlation effects, and use the simple Hartree 
approximation. We do this because of the extreme difficulty of making 
any estimate of the contribution to the elastic constants from the correlation 
energy, which in metals of valency higher than one must depend on the 
zone structure. It would in principle be possible to make an estimate of the 
contribution from the exchange energy, which must also depend on the 
zone structure, but as these two terms are so closely related to each other it 
would be difficult to interpret the significance of the results. 

In the Hartree approximation, the total energy W per atom, referred to a 
system of separated ion-cores and valence electrons, is given by the equation 


e  rp(r)p(s) 
NW=23,4,— 5 i [Free drds+ 22,U (RyRy), (Il), (1) 


where the summation over k extends over all occupied states, the integra- 
tions, and summations over | and 1’, are taken over the space occupied by 
N atoms, and the «, are the one-electron energy parameters in the Hartree 
equations 


2 
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Here V(r) is the self-consistent field given by the equation 


= 2 [_p(S) 
Vie)= $e) +e | EEE as, 
where ¢(r) is the potential energy of an electron in the field of all the ions 
in a given volume, p(r) is the number density of valence electrons given by 
the relation p(r)=22,|,(r) |2, and the integration is taken over the same 
volume. The y, are supposed normalised over the volume occupied by the 
N atoms. R, denotes the position vector of an ion, and UR) is the 
energy of interaction between two ions at a distance R apart. 

Let us now define an energy E(k) for every state by the relation 
€,—=€)+H(k), and let us introduce the wave function of the lowest state 
normalized over a single atomic cell, viz., u(r), where Up(t) = V Ny (r). 
If we then use equation (2) for k=0, and let Z denote the number of valence 
electrons per atom, we find that NW is equal to the sum of the following 
four terms : 

(i.) Ban 7 cee { UV 2p At, 

cell 


822m 
the kinetic energy of the lowest state multiplied by the number of electrons ;. 
(ii.) 22, E(k), 


the Fermi energy ; 
(ii.) — SRZ~U,(| R—Ry |) +2 | dtejagleyde 
Ze? ¢ uz(r)u2(s) ; 
+ [Arash ards, dt’, 


which is the potential energy of a system consisting of the ion-cores and a. 
charge distribution of density —Zeug(r); and 


ts [Pate coterie) “re drds, 
2 


(iv.) ir—s| 


which is the negative of the self potential energy of the charge distribution 

—e[Zu2(r)—p(r)] which ae eee in the case when the 
imation %,—y, exp(vk.r) 1s vad. 

oO ee eee the contributions made by these terms to the 

elastic constants. ; 

Calculations by Raimes (1950) of u in magnesium, using the sphere 
approximation of Wigner and Seitz, indicate that in divalent as well as 
monovalent metals the wave function of the lowest state is very flat near 
to the boundary of the atomic cell. The first stages of a similar calculation 
for aluminium indicate that the same is true in this case. It follows, as has 
been shown in the case of monovalent metals, that there is no gua 
contribution to the elastic constants from the kinetic energy of the lowest: 


state. 
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The term in (iii.) which represents the ion—-ion interaction may be 
separated into a Coulomb term, and non-Coulomb terms arising from 
polarization effects and exchange repulsion between closed shells. An 
estimate of the clearance between ion-cores, as given by Jones (1949), 
shows that in the case of aluminium contributions to the elastic constants 
from these non-Coulomb terms are negligible. 


§ 3. THE CoULOMB TERMS. 

Let us now consider the contributions to the elastic constants from the 
term (iv.) above and the Coulomb part of (iii.). With regard to the latter, 
there is nothing different from the case of a monovalent metal. We assume 
that the distribution w2(r) around each lattice point moves, effectively 
unchanged, with that lattice point as the metal is sheared. Then the 
contributions to the elastic constants are approximately the same as those 
from a set of positive charges ZeQu3(r,) in a uniform sea of negative 
charge of density —Zeu2(r,), where Q is the atomic volume in the metal 
and w2(r,) is the value of u3(r) at the boundary of an atomic cell. 

We now consider whether it is possible to estimate the contributions 
from the term (iv.) in a similar way. As the boundary conditions of those 
wave functions which correspond to states on the bounding planes of the 
Brillouin zone require that they vanish over certain planes of the atomic 
polyhedron, we should expect that the value of p(r) at a point on the 
boundary of the atomic polyhedron would be less than Zuj(r,). The 
distribution —e[Zu?(r)— p(r)], which is electrostatically neutral as a whole, 
has thus a resultant positive charge in a region about a lattice point, and 
is negative outside this region. Let us now consider what happens as the 
metal is sheared. Since both the self-consistent field and the shape of the 
atomic polyhedron change, it follows that the wave functions, and so 
p(r), will change in a very complex way. p(r) has, however, two general 
properties, apart from its periodicity, which depend only on the symmetry 
of the metal lattice; its normal gradient over the boundary of an atomic 
polyhedron is zero, and it has in the equilibrium position, in common with 
V(r) and ¥%,(r), the full cubic symmetry of the lattice. It is thus, like 
us(r), very nearly spherically symmetrical near the centre of an atomic 
cell. It is reasonable to assume that this nearly spherically symmetrical 
part moves almost unchanged with the lattice point. It also seems 
reasonable, as a first approximation, to average over the surface of the 
atomic cell and neglect actual details of what happens near to and over the 
different parts of the boundary. The contributions to the elastic constants 
are then approximately the same as those from a set of point charges 
e22[Zus(r,)—p(r.)] in a uniform sea of negative charge of density 
—e[Zus(r,)—p(7,)], where pls) is the mean value of p(r) over the boundary 
of an atomic cell. 

The contributions to the elastic constants from a face-centred cubic 
lattice of point charges +-e set in a uniform sea of negative charge of density 
—e§2 have been evaluated by Fuchs (1936), and are 


2(C,,—Cy)=0-1058(e2/2a) and C,,=0-9479(e2/2a) per atom, 
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where a is the lattice constant in A.U. From our discussion above it 


follows that we have to multiply these values by a factor Ze given. by the 
relation : 


Zeg= 2 {[Zug(r,)?—[ZaR(r,)—p(r,) P}. 
In aluminium, ug(r,) is almost equal to the average value of w2(r), so that 
putting Z=3 we have 


Z2~=9—[3— Qp(r,) bec. MR okra) 


Finally, inserting the value of the lattice constant for aluminium 

{a=4-041 A.U. at room temperature), we have for an estimate of the 

contributions to the shear constants from Coulomb terms : 
3(C,,—C,)=Z2,.0-183. 1014 dynes cm.~2, 

and C,,=Zé_. 1-639. 101 dynes cm.~, where Z2, is given by (3). 

We see at once that this contribution to C,, is already by itself much 
greater than the observed value. It is difficult to see how the value of 
Qp(r,) could be less than say 1:5, and this gives a minimum value for 
Z?.¢ of 6-75, and corresponding contributions to the elastic constants of 
3(Cj,;—C,.)=1-2.104 dynes cm.-2 and C,,=—11-1.104 dynes cm.-2, to be 
compared with the observed values of 2:1 and 2-9.1011 dynes em.~2 
respectively. 


§ 4. THe Fermi ENERGY. 


Since we have no exact knowledge of the energy surfaces in momentum 
space, and moreover no idea of how they change as the metal is sheared, 
we can only proceed by making simple assumptions consistent with known 
general properties of these surfaces. 

Aluminium has a face-centred cubic structure and therefore a body- 
centred reciprocal lattice, and the Brillouin zone is the well known 
truncated octahedron. The energy surfaces cut the bounding planes of 
the zone normally ; they are also generally considered to be approximately 
spherical near to the origin, and spheroidal in the second zone near to the 
centres of the faces. “There are three valency electrons per atom, the first 
zone containing sufficient states for exactly two of these with each state 
doubly occupied. The remaining electrons occupy states overlapping the 
hexagonal faces, which are the nearest- to the origin, with possibly a few 
overlapping the square faces. There is also the possibility that the first 
zone may be not quite full. 

Let us consider how the reciprocal lattice, and so the Brillouin zone, 
changes as the metal undergoes a shear. The strain energy W of a cubic 
lattice is given in terms of the three independent elastic constants, C,,, 


C,, and C,,, by the equation 
W=W,+3C, 1 (e2be2+ 62) +Cyo( Cyl + eee tently) + 2CaalYyet Yao t Vays 
. Beet) 
where ¢,, ¢, and e, are the diagonal components of the strain tensor €, and 
$V ye BY x2 aNd 47 are the shear components. We shail choose two pure 
shears which correspond respectively to 3(C,;—12) and C,,. 
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A shear corresponding to 4(C,,—C,2). Consider a strain consisting of 
compression or extension in the (001) direction accompanied by that 
uniform expansion or contraction in the perpendicular plane necessary to 
keep the volume constant. We can express the sheared face-centred 
cubic lattice as follows : 


al 1 a l 1 al 
— j —_ = - =- —— —_ i 1, 0), 
a= 5 a(” iF 2) ao Sa(ho 2) a3 5A | ? ) 


where in the equilibrium position A= é=1, and to keep the volume constant 
during a finite strain we must take A?=1/€. “Remembering that the 
displacement t at r is given by t=r.e, it follows that the components of € 
in this case are as follows : 


e,=ey=1/A—1, e—=1/€A—1, and y.—Ya=Yay—0- 
Inserting these values into (4) we find that 


7 ; 3 (d* is 
MCu—Ci)= iG ut eS ee 
where the suffix zero implies the value at the equilibrium position. The 
sheared reciprocal lattice, defined by the relations a;.b;—276,;,, where 5;; 
is the Kronecker delta, is given by the vectors b, as follows : 


In Qn In 
b=—A(-1,1,8, b=—A(,—-1,8, b= ACL, 1, —€). 


yw 


A shear corresponding to C,,. Consider a strain which is the same as the 
previous one except that the (001) is replaced by the (111) direction. The 
metal lattice in this case may be expressed in the form : 


aa aé 
Ais Bl. Lb 1), a= 5 All. a—l, 1), as= 5 BRU 1, y—1), 


where 6=2[(1+-)(2—»)]-?, and in the equilibrium position B=8=y=1. 
For constant volume we must take B°=4[(1+-7)?(2—n)]4. In this case 
we have ¢,=¢,=e,=}3(8/B)(3—y)—1, yy=Y2=Yay= (8/B)(n— 1), and 


C 1 /d?W 
u=3(Gr),: e. « «e «= 36 ae (6) 


The reciprocal lattice is given by the vectors 
__ 2n 2a 2Qar 
bh B (—», 1, 1), a eae B (1, —9, 1), b—— Ba, 1, ap] 


In each case of course the first derivatives of the energy vanish at the 
equilibrium position. 

The faces of the Brillouin zone are formed by the planes which bisect 
perpendicularly the lines joining the origin of the reciprocal lattice to its 
nearest and next nearest neighbours. We can thus determine in each of 


the two strains the altered positions of these planes and of the corner 
points where they meet. 
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We shall consider separately the contributions to the elastic constants 
from a completely filled first zone and from overlapping electrons. We shall 


also consider what the effect of holes in the first zone would be if th 


ey were 
present. 


Contributions from a completely filled first zone. The following treatment 
is suggested by regarding the actual energy surfaces as the spherical 
surfaces given by the relation 

h?k? ‘i 
k) => Xo SE . ° ° ° ° F ° . ( ‘ ) 
so modified as to cut the faces of the zone normally. We find that when 
we approximate to these modifications in a simple empirical way the 
resulting contributions to the elastic constants differ only slightly from the 
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A typical tetrahedron, in this case at a square face. 


values we would obtain by considering simply the change in the energy of a 

zone in which the surfaces are given by (7). We ae that this is true 
to the faces. 

whatever the exact nature of the surfaces near ; 

Let us suppose that the first zone energies at the centres of all the faces 
and at the corner points are the same fraction (1- -A) of the energies at these 
points given by (7). Let us also suppose that % and A remain die 
during a shear, so that the energies at these pony as as ee 

i igi then divide the zone into tetrahe 
the distances from the origin. Let us 
(fig. 1) whose vertices are the origin O, the symmetrical centre RP iif ik 
(which even in the strained zone is the point at which it is we <a : 
perpendicular from the origin), a corner point R “a ae eee anc ‘ ie 
i f the face adjacent to R. 
endicular from P to an edge o Bes 
= ey of these tetrahedra, 12 to each hexagonal face and 8 to 
each square face. 
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We have now to find an empirical relation which expresses E(k) as a 
function of k within a single tetrahedron and satisfies the following 
conditions. If taken in all the tetrahedra with the appropriate dimensions 
in each one, it must fill the whole zone with energy surfaces which are 
continuous and continuously differentiable. We require it to satisfy the 
assumptions we have made, and it must possess the general features we 
have already described. Such a relation is 


h2 ki. 2/A k 2/A k 2/A 
rarnatic(e[o) oe) eT 


where k,, k, and k, are measured in each tetrahedron in the directions of 
OP, PQ and QR respectively, and p, g and r are the appropriate lengths of 
these edges. Fig. 2 shows how the parameter A modifies the surfaces so 
that they cut the faces normally. E(k)/(«,h?—p?/87?m) is shown as a function 
of k,/p along OP for different values of A. 


Hig..2: 

1-0 A=0 
A=0: 
A=02 
A=03 

0:5 
0 —_| 
0 0:5 1-0 


‘The effect of the parameter A in expression (8). K(k)/(agh?:p?/872m) is shown 
as a function of k{/p along OP. 
We have now to use these surfaces to calculate the first zone part Wi. of 
the Fermi energy. The contribution from a single tetrahedron is 
2 I 2 5 yest) 
373 “ga Parl Fp* + Gq?+Hr*], 


‘where 
hee a ag ds 
10 2(2+3A)’ ~~ 20s 2025 -SA\(IEEA) * 
1 MM 


~ 60 22FBNT-FANA)’ 
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For each tetrahedron p, g, and r are functions of the strain parameter 
€or. We find that because of the high symmetry of the distortions ee 
have chosen there are only six different tetrahedra in the first shear, and 
five in the second. Summing over all the tetrahedra we obtain wr as a 
function of the strain parameter in each case. For the shear Op Sesrne 
to 3(C,,—C,,) we find for the contribution to the energy per unit volume 


h2 
ss $9 F(16€ + 365— 165+ 3é7) + G(16é— 1248+ 16£5— 5é7) 


wi 
+H(16g—20€9+ 8£54 27)], 


and for the shear corresponding to C,, 
3a ,h? 4 a | : 
es La ee CANS 5 
F~ 9ma® Resicel (1-49) [F(16 1287-4 406n?—341y 
+ 96n*-- 109° — 6°—n") 
G(40 3007+ 5787? — 49973 + 18474—22n5—2n6+-7) 
+ H(64 —152n+ 15872—1217?-+ 96y!—46y°+ 107°—n?)]. 


From these expressions we find that 


(Fe) = (=) =i) 
dé Sie dn on. : 


-and, using relations (5) and (6), that the contributions to the elastic 
constants from a completely filled first zone are 


3(C,3—Cyy) = ae 75 ia Gb WV Ff PLD WY 3 nl WY 


ah? E 6)2(4-+2A—1122) ] 
3ma> (2+3A)(1+A)(2+A) |’ 
-and 
gh? 28-72] 
C8 ina’ [ © (@+3A)(1-FA)(2+)) |’ 


We see that the effect of taking a value of A different from zero is very 
-small. Even with a value of \ as large as 0-3 the variation in the elastic 
-constants is only of the order of 0-3 x 1014 dynes cm.~*. We shall therefore 

take as our estimates of the contributions to the elastic constants from a 
filled first zone the values corresponding to A=0. These are 


1(C,,—Cy2)=1-5%. 10% dynes em. and Cyy=4-5a9. 10% dynes em.~ 


-at room temperature. We do not know the value of x) in aluminium, but 
we are unlikely to be involved in serious error by taking it to be unity. 


Contributions from overlapping electrons. We shail make here the 
_simple assumption that during shear the energy surfaces over a face in the 
second zone move with that face without change relative to it. We shall 


_also neglect the small effect of finite temperature on the Fermi distribution, 


and suppose that the states are completely filled with two electrons each 


‘up to an energy ¢. 
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It is easy to see what happens as the metal is sheared, if we remember 


that the number of electrons per atom remains exactly three, and the 
energy ¢, while changing its value, is the same for different parts of the 
Fermi surface. As the distances of the faces from the origin vary, electrons 
will transfer from those faces receding from the origin to those approaching 
it. The change in the Fermi energy has thus two causes : the shift of the 
energy surfaces as a whole with the faces, and the relaxation effect of the 
electron transfer. Let us examine this in detail. 

It follows from the symmetry of the reciprocal lattice that the perpen- 
dicular distances from the origin of opposite faces of the zone are always 
equal, so that we may take these faces together in pairs. Since 
faces of one kind, hexagonal or square, do not in general remain at 
equal distances from the origin during shear, we shall refer to the pairs of 
faces by two suffices, 7 denoting the kind of face, and j the different pairs 
of that kind. “Let us define the density of states N(E) as the number of 
electrons per atom per unit energy range, allowing for two electrons being 
in each state, and denote by N,(E—E,,) the contribution to N(E) from the 
pair of faces with energy E,,; at the centre of a face. Then the total number 
of electrons per atom, n,;, occupying states over this pair of faces in the 
sheared state is given by 


¢—E; 
ng= | uN (c)de xe ne 


If we denote by Wi! the total contribution to the Fermi energy from 
overlapping electrons, we have that 


Wr=2; ; Wi: 
where We= |. N,(E—E,,)E dE 


2] 
Wy 


c—Ei; 
=nyEyt | N.(cje de. 2) or 


If we let a stand for either strain parameter, € or y, and differentiate 
equation (9) with respect to x, we have 


dns dt dk, 
Te = Nul—Ey) ( — a): I i i 
Since the number of electrons per atom does not depend on x, we have 
dn; 
Pi perc ith > Chane EMME, (00) 


If we assume that during shear the energies at the centres of all faces of one 
kind, hexagonal or square, are the same function of the perpendicular 
distance of the face from the origin, it follows that dK,;/dx is proportional 
to dp,;/dx, where the constant of proportionality depends only on 7, and so 
2',dK;,/dx is proportional to X;dp,,/dx. We-find below that the vanishing of 
2ij(Upjjld2)g for all 7 follows as a geometrical consequence of the volume 
remaining constant. We have then 


dE,, 
2; (=), ——()5 . . . . . . . . (13) 
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‘Taking together equations (11) and (12) in the equilibrium position, and 
using equation (13), we find that 


d¢ 
(i) =0. RPA eat sin 


If we now differentiate equation (10) twice with respect to a, use equations 
(11), (12), (13) and (14), and sum over i and j, we find that (dWH/dzx),=0, 


and 
ewu PE, JE.\2 


0 
where n;, Cy and E, are the equilibrium position values of n;, ¢ and H,, 
respectively. 

We can obtain estimates of the values of the summations over yan 
equation (15) by assuming that the energies E,; vary as the squares of the 
distances from the origin of the points they refer to. Equation (15) then 
has the form 


wi - 
(SS), =2,[A.nE;—yNj(Co—E,) EH], . . . . (16) 
where the A; and ,; are numerical coefficients. 

The squared distances of the faces from the origin in each shear are as 
follows. In the shear corresponding to 4$(C,,—C,,), the four pairs of 
hexagonal faces are at the same distance (z/a)A(2+ &)# from the origin, the 
pair of square faces in the (001) direction are at a distance (27/a)A€, and 
the other two pairs are at a distance (27/a)A. In the shear corresponding 
to C44, the pair of hexagonal faces in the (111) direction are at a distance 
(z/a),\/3B(2—n) from the origin, the other three pairs are at a distance 
(z/a)B(2+-7?)*, and the three pairs of square faces are at a distance 
(7/a)1/2B(3—2n+7?)?. From a knowledge of these distances it is easy to 
evaluate the coefficients A, and p; for each shear, and to verify the result 
»,(dp,;/dx))=0 in each case. The following table shows the contributions 
to the elastic constants per atom from overlapping electrons over 
hexagonal and square faces, obtained from equation (16) using relations 


(5) and (6). 


Hexagonal faces Square faces 
$(Cy,;—Cy2) Sg Hg 2np, Hy —2N, kz 


re aes) RS Ae 5) eee 
Nq and n, are the numbers of electrons per atom in the equilibrium ane 
over one pair of hexagonal and square faces respectively, N, and ‘i the 
corresponding contributions to the density of states per atom at the see 
surface, and E, and E, the equilibrium position energies at the oe es 6 
the faces. In each case the positive term involving the A coefficient iy 
be regarded as arising directly from the shift of the energy surfaces as a 
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whole with the faces, and the negative term involving the p coefficient as 
arising from electron transfer between the faces. This electron transfer 
term only arises of course when some of the faces in question approach the 
origin and others recede from it. Thus it does not appear in the contribu- 
tions from hexagonal faces to }(C,,—C,,) or from square faces to Cy,, as 
in these cases the faces remain at an equal distance from the origin. 

We can further reduce the number of parameters involved by approxi- 
mating to the energy surfaces with spheroids. For the surfaces over the 
face given by the sufficies i and j, we have 

2 
BABy + ge (EET KE toy ky) 

where k, is measured perpendicular to the face, and «; remains constant 
during shear. We then have 

Q 47 1 [82m 3/2 
HBF aE]. 

3. ON; a 
and Nees = te MPR ee ee 

The effect of the presence of holes in lhe first zone. We can proceed 
exactly as with overlapping electrons, the only differences being changes of 
sion. The following table shows contributions to the elastic constants 
corresponding to two ways in which holes could occur, viz. (i.) at the corner 
points, and (ii.) in approximately plane regions over the square faces, as 
suggested by Matyas (1948). 


nN 


Holes at the corner points | Holes over the square faces |, 


2N_E2 


% 6.651 o5-,6uc 


n, is the number of holes per atom in one closed region, formed by the holes 
at four corner points (there are six closed regions in all), and n, is the 
number over one pair of square faces; N, and N, are the corresponding 
contributions to the density of states per atom at the Fermi surface; and 
E, and EK, are the first zone energies at the corner points and the centres of 
the square faces respectively. It will be noticed that the terms arising 
from the shift of the energy surfaces are negative for holes. 


§ 5. Discusston. 


As we do not know the values of any of the parameters which appear in 
our results, viz. p(r,), «, and «,, E, and K, etc., we cannot make a direct 
calculation of the elastic constants. We can however see whether it is 
possible, by assuming reasonable values for these parameters, to account 
for the observed values of the constants as the sum of the three contribu- 
tions which enter into the present theory. In doing this we can make use of 


a 


ee EEE 
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two observed quantities which depend on the band structure. These are: 
(i.) the width of the soft X-ray emission spectrum, given by Skinner (1940) 
as 11-8 eV., which corresponds to ¢,; and (ii.) the electronic specific heat, 
given by Kok and Keesom (1937) as 3-484 x 10-4T cal. deg.—mol.—, which is 
given in terms of the density of states per atom at the Fermi surface by the 
expression (7°/3)kKN(Cy) RT, giving N(¢,)=0-614 per eV. per atom. 

Let us leave aside for the moment the possibility of holes in the first 
zone. If there is only overlap over the hexagonal faces, then we have 
4n,—=1 and 4N,—N({,), and the contributions from overlapping electrons 
to 3(C,,;—C;,2) and Cy, are 3E, and 3E,—4N(¢,)E? respectively. Using 
relation (17) we find (,—E,—2-4 eV., and therefore E,—9:4eV. These 
contributions to the elastic constants are then }(C,,—C,,)=6-3 eV. per: 
atom=6-1 x 104 dynes cm.~?, and C,,=— 17-8 eV. per atom=— 17:3 x 1014 
dynes cm.-*. The very large negative contribution to C,,, due to electron 
transfer between the hexagonal faces, is just what we are looking for to. 
balance the large positive Coulomb and first zone contributions. The: 
contribution to $(C,,—C,,) however is even by itself much greater than 
the observed value of this constant, viz., 2:1 x1044dynes cm.. This. 
value can only be reduced by supposing that there is some overlap over the 
square faces as well as the hexagonal faces. We then have a negative 
contribution to $(C,,—C,,) from electron transfer between them. We now 
have five parameters, viz. p(r,), H,, E,, ~,and«,, to relate to four observed. 
quantities, viz. f,, N(¢,), and the two elastic constants. Since the number 
of overlapping electrons over the square faces is small, it is reasonable to. 
use the approximate expression for «, derived from the theory of nearly 
free electrons, viz. «,—1+2E/(E,—E), where E{” is the free 
electron energy at the centre of a square face, ¢. e. (h?/87?m)(27/a)?=9-2 eV. 
We then find that the values we have to assign to the parameters, in order 
to reproduce the observed values of the elastic constants, are as follows : 
Z?,—=6-9,, and therefore Qp(r,)=1-5g; E,=9-0,eV. and a,=1-9,; H,=11-2, 
eV. and «,=10-2,. The corresponding values of n, and 7, are 0-24, and. 
0-01, per atom respectively, and those of N, and N, are 0-133, per eV. per 
atom and 0-027, per eV. per atom respectively. The following table. 
shows the constitution of the observed elastic constants as explained by 
this model. The units are 104 dynes cm.~?. 


Contributions 4(Cii-C,3) Ca | 
| 
Coulomb 1-3 11-4 
First zone 1-5 4-5 
Overlapping electrons —0-7 —13-0, 
Total (observed) 2-1 2-8, 


It will be noticed that the value we have to assume for «, is much lower 
than would be given by an expression corresponding to the one we have 
used for «,. This is not a serious objection, however, as the region of 
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occupied states over the hexagonal faces goes well beyond the range of 
validity of such an expression. On the other hand it still seems reasonable 
to approximate, as we have done, to the general shape of the energy 
surfaces with spheroids, while leaving «,, unspecified. 

It is of interest to enquire whether the observed constants could be 
explained by assuming the existence of holes in the first zone instead of 
overlap over the square faces. We see however that in the case of holes at 
the corner points the coefficient of the term in $(C,,—C,,) arising from 
electron transfer is very small. In fact, if we assume spherical energy 
surfaces about the corner points, we find that in order to account for the 
-observed values of both elastic constants consistently with known values of 
¢, and N(¢,), and with a value of Qp(r,) not less than 1-5, it is necessary to 
_assume that there are at least 0-3 holes per atom (7. e. 0-05 in a single closed 
region), and that the corner point energy E, is at least as high as 16-5 eV. 
This is 5 eV. above the free-electron energy of (h?/87?m) 5(7/a)?=11-5 eV. 
In the case of holes over the square faces, the coefficients of the 
terms are of course the same as for overlapping electrons there, the 
-only difference being the sign of the term arising directly from the shift 
of the energy surfaces, and this is very small as long as the number of 
-electrons or holes is small. It is thus possible to explain the elastic 
constants in the same way. We should however have to assume a value 
of E, a little greater than ¢), 7. e. about 3eV. above the free-electron 
energy of 9:2eV. These high first zone energies would imply energy 
surfaces which are very different from those of the simple model we have 
been using, and the present theory is not really adequate for their discussion. 
It is clear, however, that the elastic constants cannot be explained by 
assuming a small number of holes at the corner points in the same way as 
they can be by assuming a small number of overlapping electrons over the 
square faces. This is because in the shear corresponding to $(C,;—C,.) 
the corner points do not move as much as the centres of the square faces, 
so that the resulting amount of electron transfer is smaller, and the 
corresponding negative contribution to the elastic constant is not large 
enough to balance the large positive contribution from the hexagonal 
faces. There still remains the possibility that there are a small number of 
holes at the corner points, as well as overlapping electrons over the square 
faces, but this would have only a very small effect on the elastic constants. 


§ 6. VARIATION OF THE ELAstic ConSTANTS WITH ELECTRON 
CONCENTRATION. 


When metals of valency higher or lower than three are dissolved 
substitutionally in aluminium, and the valency electrons of the soluble 
atoms are shared with the whole metal, there is a change in the number of 
occupied states in the band scheme. The consequent change in the height 
of the filled region leads to a change in the contributions to the elastic 
constants from overlapping electrons. The full curves in fig. 3 show the 
effect on the total elastic constants of the variation of these contributions 
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with the electron concentration per atom. The large rise in $(C,,;—C,.) 
with falling electron concentration is due to the decrease in magnitude of 
the negative term arising from electron transfer between the square faces 
in the corresponding shear. This term is proportional to the contribution 
from these faces to the density of states at the Fermi surface, and so is 
proportional to the cube root of the number of overlapping electrons there. 
There is thus a very sharp change in $(C,,—C,,) at an electron concentra- 
tion of approximately 2-67 per atom, which is the point at which there is 
just no overlap over the square faces. Similarly, as the electron concen- 
tration falls still further, and the number of overlapping electrons over the 
hexagonal faces becomes small, C,, begins to rise rapidly. There are also 
changes in the details of the energy surfaces in momentum space, and so in 


Fig. 3. 


Variation of the elastic constants with electron concentration (units 
1011 dynes cm.~?). The full curves show the effect of the variation of the 
contributions from overlapping electrons alone. In the broken curves 
an allowance has been made for the change in the Coulomb contribution. 


the density of states curve, but we should expect the effect of these changes 
on the elastic constants to be relatively small. 

At some point as the electron concentration decreases, holes will appear 
in the first zone. From the results we have obtained in §4 for the 
contributions to the elastic constants from holes, we should expect that 
this would have little effect on C,, but would lead to a sharp fall in the 
value of 4(C,,—Cj,). The magnitude of this would depend on the 
curvature of the energy surfaces near the point of highest energy in the 
first zone. If there were no overlap over the square faces in pure aluminium, 
but instead holes in the first zone, then the variation of the elastic constants 
with electron concentration would be different from that shown in fig. =F, 
and would depend on the position of the point of highest energy in the 
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first zone. If the holes were at the corner points there would be a rather 
large decrease in both }(Cy,—C,,) and Cy. with decreasing electron 
concentration. If they were over the square faces the effect for small 
changes of electron concentration would be similar to that shown in fig. 3; 
it is not possible to say what would happen at lower electron concentrations, 
except that we should not expect a sharp change in 3(C,;,—C,») of the 
kind that arises from the disappearance of overlap over the square faces. 


Fig. 4. 
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The contributions to the density of states, N(E) per eV. per atom from 
overlapping electrons over (a) hexagonal faces, and (b) square faces. 
The full curve is the sum of curves (a) and (6). 


As well as changes in the Fermi energy contributions to the elastic 
constants, there are also changes in the Coulomb contributions arising 
from the changes in electron density in the metal. This effect is probably 
small in the case of $(C,,—C,,) because of the small Coulomb contribution, 
but it is important in the case of C,,. In order to obtain an estimate of the 
size of this effect in the case, for example, of a divalent solute, let us 
consider a crystal of the pure solute in a hypothetical face-centred cubic 
phase with the same lattice constant as aluminium. We should expect the 
value of Qp(r,) in this crystal to lie between say 1-0 and 2-0, and so Z2 
between 3-0 and 4:0. As the concentration of this solute in cient 
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varies from zero to unity, the Coulomb contributions to the elastic 
constants change from the values for pure aluminium to those for this 
hypothetical crystal. The values between these two limits will depend on 
the exact nature of the changes in the electron density. The broken 
curves in fig. 3 have been obtained by supposing that Z2, varies linearly 
with solute concentration between the values of 6-9, Or pile aluminium: 
and 3-5 for an electron concentration of 2-0 per ont. It is difficult to see 
how the net effect on the elastic constants can be such as to give rise to 
values which are very different from these. 

There does not appear to be any experimental data available relating to. 
the change in the single crystal elastic constants with electron concentra- 
tion. One difficulty is that only very small quantities of other metals are: 
soluble in aluminium at room temperature without a phase change occurring. 
The most soluble are magnesium and zinc, which dissolve to the extent of 
about 2 atomic per cent, corresponding to a fall in electron concentration of 
0-02 peratom. This should give, according to the broken curves in fig. 3, 
an increase in $(C,,;—C,,) of approximately 0-05.1014 dynes cm. and a. 
decrease in C,, of approximately 0-12.10! dynes cm.~*. This corresponds 
to a change in the ratio C,,/4(C,,—C,.) from 1-37 to 1-28, a decrease of 
6-4 per cent. The only way to measure the effect of larger changes of 
electron concentration would be to make observations either at higher 
temperatures or on quenched solutions. The most favourable case for 
this is that of zinc, which has a maximum solubility of 66-5 atomic per 
cent at 382°C. This corresponds to an electron concentration of 2-335 per 
atom. Other possibilities are magnesium, with a maximum solubility of 
16 atomic per cent at 450°C., and silver, with a maximum solubility of 
18 atomic per cent at 538°C. These correspond to electron concentrations 
of 2-84 and 2-64 per atom respectively. Copper has a maximum solubility 
of 2-5 atomic per cent at 548°C, corresponding to an electron concentration 
of 2-95 per atom. Silicon appears to be the only element with more than. 
three valency electrons per atom which is appreciably soluble even at. 
high temperatures, the maximum solubility in this case being approximately 
2 atomic per cent at 577°C. This corresponds to an electron concentration: 
of 3-02 per atom. 
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ABSTRACT. 


Schubert (1946) and others have drawn attention to the fact that the 
‘standard derivations of the formule for the behaviour of a perfect Bose- 
Einstein gas are invalid in the region where condensation should occur. 
In this paper a new and rigorous derivation is given of formule for the 
macroscopic behaviour of such a gas. 


§1. INTRODUCTION. 


ScuuBert (1946) and Dingle (1949) have pointed out that discussions 
of the phenomenon of condensation in a perfect Bose-Einstein gas by 
London (1938) and Fowler and Jones (1938) are invalid because of the 
divergence of some of the series used. Temperley (1949) has considered 
the behaviour of a model with equally spaced energy levels and has 
concluded that the usually accepted formule are incorrect in this case. 

In this paper formule for the macroscopic behaviour of a perfect 
Bose—Hinstein gas will be found by a method similar to one used by Kahn 
and Uhlenbeck (1938) in considering the condensation of an imperfect 
gas. We shall have in mind the particular case of a perfect gas in a cubic 
box, but the results will be expressed in a more general form so that they 
may be applied to a perturbed system such as a perfect gas in an external 
field. 

In a second paper (Fraser 1951) formule for the occupation numbers 
of the energy levels will be obtained. 

All the thermodynamic functions for a perfect gas of N particles can 
be found from the partition functiont 


Ly=d" exp {peng}, . . ° . Si é ‘ . (1) 

v 
where .=1/kT, 2” means the sum over all values of the set of positive 
integers %, 21, NM, . . . such that 2n,—N, and ¢, «,, €,, . . . are the 


vu 
energy levels for the individual particles. In particular, the free energy 
of the gas is given by the relation 


PS WT log Ze eee le (es. 


* Communicated by the Author. 
+ See, for instance, Schroedinger, E. Statistical Thermodynamics (Cambridge 
1945). : 
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and all other thermodynamic functions can be derived from this by well 
known relations. 


it a Bose—Hinstein gas it is easily seen that Z, is the coefficient of 
Grin 


ey lee xP (—pe))-+0 exp (—2ye) +... P= —€ exp (—He), 
(3) 


where w; is the weight of the level with energy «; We can therefore 
write 


dos = sa hC-OPTTL —L exp (—pa)] de, 9 surety 


where the contour includes the origin, but excludes the poles C= exp (ie;): 
on the real axis. In the method used by Fowler (1936) this integral is 
evaluated by the method of steepest descents, but this has been shown 
by Schubert to be invalid in the neighbourhood of the condensation point. 

In this paper we shall find an expression for the free energy per particle 
as the volume and the number of particles tend to infinity in such a way 
that the number per unit volume remains constant. 


§2. THE ASYMPTOTIC ForM OF THE PARTITION FUNCTION. 


The energy levels «; will depend on the volume of the gas and we must 
first impose some restrictions on this dependence. We write 


= Zepexp(—sue)=a,(V), bie EM LN) ety! (5)i 


where a,(V) is a function of the volume only, and require that there 
exists a set of functions b,(p), depending only on the density p=N/V, 
such that, for N>M, where M is any integer, 


(i) a,(V)<(1+8y)NO,(p) for all s, 
(ii) a,(V) >(1—45y)Nb,(p) for all s<s8y, 


(iii) the radius of convergence of 2 b,(p)a* is unity, 
s=1 


and (iv) 2'sb,(p) is finite, 
. s=l 


where 5, tends to zero and sy tends to infinity as M tends to infinity. 
We will in future write merely b, for b,(p). A discussion of these: 
restrictions will be relegated to the appendix where it will be shown that 
they are satisfied by the energy levels for a perfect gas in a cubic box.. 
We see from the definitions that the b,’s are always positive. 

Rewriting (3), we find that Zy is the coefficient of gX in exp 2a,C", 


r 
and hence that Zy may be considered as a monotonically increasing: 
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function of the a,’s. It follows that if we define functions Zx y and 
Ly, y Such that 


1 foe) 
Zhu = ag GOO exp {N+ Oy) 2 ddtp ad, 


6 = hE exp IN(I—3y) 2°b,6° hdl ) 
Zyu= 56 exp 4 N(1— mu) & s ; 


then Zy y>Zy>Zy, » for all N>M. 


The contours are to be taken round the origin but inside the radius of 
convergence of = b,¢%. Zy 4 and Zy y are always positive since the 
5's are. 

We can now find 


1 \ 
(es ats /, / hee wks PS he 
log R’= ls (5 log Zy, x) and log R ut (5 log Zy, x) 
by observing, following Kahn and Melons (1938), that R’ and R” are 
the radii of convergence of the series YZ y,a% and YZy yx%. In these two 


series, Zy y and Zx y refer to the same density for all N, while in generating 
series of the type (3) the Zy refer to a fixed set of energy levels and there 
fore to a fixed volume 

We can now transform the functions YZy yeX and 2Zy ye% in a 
manner which is precisely the same as that of Kahn and Uhlenbeck 
(1938) and which we shall consequently only outline here. 

Using (7) we may write 


SZy yx® = sa TOO EOD exp {N (1+8y) Zb bf ‘hat 
N 7d Ad 


= mae pees Sete “a ae(S) 
{2 exp {(1-+54) é vcs} 


so long as we choose the contour so that 


exp {(1-+5y) ) Zo bcs} 
= 


everywhere on the contour. This condition is sufficient for the validity 


of the interchange of the sum and integral operations, and can always be 
satisfied for small enough «. 


For sufficiently small « it may be shown that 


ms <l—e, e>0, 


(x exp {(1+y Fo ah 
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has only one root, Cy, say, inside the contour. This is the only pole 
which contributes to the integral (8). We can therefore write 


1 1 
SEG i) Dole WRAY = fee 
Ss = hse epa) eae ®) 
where fre) =p exp {(1-+8y) 25,0}, oe GS LUO) 


1 exp {(1-+3y ) Bb, af 
= 5 rarer aeons . . . ° . (12) 
d at, 
and <= 2 93 gal) | eee eee eg = C13) 
M 


Since we have shown that (9) is true for small a, it follows that, if the 
series in (10), (11) and (12) are replaced by their NS continuations, 
(9) will be the analytical continuation of 2Zy yx. 

Since all the Zy y are positive, the radius ot convergence R’ is such 
that «=R’ is the first singularity of 2 Zy yw as x increases along the 
positive real axis. 

XZ, yx and its analytical continuation will be regular in that region 
of the x-plane which is mapped by the region of the & -plane in which 


(2) 91, u(Su) Al 
and (b) gy w(x) is regular. 


(14) 


Conversely, at points corresponding to those at which (a) or () is violated, 
2Zy, yt™ and its analytical continuation will have a singularity. 

We see from (13) that &, is a monotonically increasing function of 
Z as x moves along the positive real axis until (a) is violated. 
Hence x=R’ corresponds to the smallest real positive value of Cy for 
which (a) or (b) is violated. Because of postulate (6, iii), the first 
singularity of g, y(¢) is at C=1, and it follows that (a) will be violated 
before (5) if and only if g, y(¢)=1 has a real positive root less than 1. 

Hence, if (1 +3y) 2186 5 ¢8=1 has a real positive root 0 y<l, 


sa (+50 2 E bi u} 
Aa 


R’ £6, M 


while, if no such root exists, 


1 — exp {aro 2b}. 
s=1 


| 
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Since Lt (Zy,y)"%=1/R’ and Zy<Zy, y for all M, we can write 


Noo 


1 1s , 
Lt (r log zx) < Lt It (y log Zy, x) = 2'b,ff—log%, . (154) 
N>o N M>o N>o N s=1 
when 2'sb,¢*=1 has a real positive root (<1, and 


1 1 ; ae 
Lt (x ee Zs) < Lt Lt (108 2%, )= 20, - + « (156) 


No M>a N>o 
if no such root exists. 
In a precisely similar manner we find that 


1 
PY ee 
aur 1—go (Cyr) 
bas 
where Go, u(S)=(1—8y) 2' sb ,c* 
s=1 
la 8u 
exp 4(1—8,,) 2 bats} 
: 1 { s=1 
and = == =>... 
% M 


In this case, since gy (2) is a finite series, YZx yx can only become. 
singular if 7, y(4.)=1. Hence 


1 7 : "sg fe 
Lt (5 log Zy, x) = (1—8y) 2 bo m— log fo, wp 
N>o s=1 
“6 (16) 
where (1—8y) 2 8b oo = I. 
=1 


Since Zy>Zy, y for all M, 
1 
Lt (5 log zx) = plat Le & log Zx, x) : 


N>o M—>o N—-o 


fo a 
and therefore we can deduce from equation (16) that, if 2 sb Ciel aiias 


a real positive root €)<1, ond 
1 L 
Lt: (= log Zy |S Bb op Er ite Slo a 
’ N>o N s=1_ 


If, however, yy sb €’=1 has no solution less than one, the solution of 


su s=1 
(1—6y) 2 sb.¢*=1 will tend to one from above as M tends to infinity, 
s=1 


since 2' sb,¢* diverges for ¢>1. In order to find the value of the limit 
s=1 
of the right-hand side of (16) as M tends to infinity, we must evaluate 


su 
Lt {0—By) 2 b65tu 
fer] 


M—o 


ony 
where (1—8y) 2 8b Co’ y= 1. 
sm 
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We write (5 w=1+ey. Since Y sb,f* diverges for £1, €y tends to 
s=1 
zero through positive values as M tends to infinity. Now 


Syp 


ee a 


$= 


su 
(1 — dy) ae b.lou= (1 —6y,) 
s= 


845 Sy 
=(1-8y){ 2 b,+ey E sb (1-+Bey)} lO =O) 
s=1 s=1 


sur 


<< (I—0yy) 2:6, eyy. 
s=1 


We also have the inequality 


Sy Su 
(One O sCorar-> (1 — dye) 
s=1 $= 


uM sy a) 
Therefore Lt {(1—x) DE bit} = Lt {(—3x) Bb,h— Dy 
M+ oo Ces M—->o s=1 61 
Hence, if 2 sb,¢*—1 has no solution less than one, equation (16) yields, in 

s=1 
the limit of M tending to infinity, 
R 1 a 
Lt (x lor x) > 20, a tosh wettest (LS) 


No N 


§ 3. ASyMpTOTIC Form oF THE FREE ENERGY. 


By combining the inequalities (15), (17) and (18) we obtain an equality 
which, combined with equation (2), may be written 
F 


= Lt —~=—AhT Lt 
f Ro Le ( 


= log Zx)=—KUY 2 b,¢i— log to} ne (19a) 
s=1 


if z sb,¢®=1 has a real positive solution Co<1; if, on the other hand, 
s=1 
there is no such solution, 
f=—kT Lt at log Zy } =—kT 26,, ween LO 0) 
Nea = 


No 


where f is the free energy per particle. 
The change in the form of the free energy which occurs when ¢)=1 
corresponds to “ condensation’. We shall show, in the appendix, that 


for a perfect gas in a cubic box 
(2armkT)3/2 
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Therefore, for such a gas, the free energy per particle is 


J a (HT)" By ea 3 TAT log fo, (T>T) 
p 
3/2 
—— ee (ET) s r ; (T<T,) 
p 
2amkT)?? y af 
‘where ore 338 =) 
and T, is defined by the equation : 
3/2 
(2amkT 9) x rae 7 


ph $=1 38/2 


These expressions are equivalent to those derived by London (1938) 
and Fowler and Jones (1938). and all the macroscopic characteristics of 
the condensation phenomenon can be derived from them. 

The expressions (19 a) and (19 b) can be written in a more familiar way 
by observing that 


1 
sb,= Lt in me exp (—spe)| ’ 


N>o 


the limit to be taken with the density kept constant. If all the series 
are uniformly convergent, we can write 


1 = (So exp (—pe;)*} 
z hair Lt tx are ane 


N->o s=1 


= Lt tr re [1—{, exp (—nea}. 


N+ 
Similar] b= Lt ae bad 
. J Zo, 3) N->o a hs ‘exp (ue;)— 1) ° 


Thus, we can write the free energy per particle, in the limit as N tends to 
infinity, in the form 


de 
f=— Ut Baeza log [1—{y exp (—pe;,)]—N log cot | , 


N-—> «© 


1 W; 


where Lt? = 3 
n->oN ; O51 exp (ue)—1 


or, if there is no solution f) to this equation, j=1. This expression 
for the free energy is identical with that which is obtained by the method 
of steepest descents. 
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APPENDIX. 


It remains to discuss the validity of the postulates (6). 

In order that postulates (i.) and (ii.) should be justified, were quire that 
the number of energy levels for each particle between E and E--dE should 
be of the form V{g(E)-+-A(E, V)} dE where h(E, V) tends to zero as V tends 
to infinity. An assumption of this nature is necessary if the free energy 
per particle is to tend to a definite value as V tends to infinity. If the 
density of energy levels has the form assumed, then 


Seo, exp (—sue)=V | {g(E)-+-h(E, V)} exp (—syE) dE 
i 


“Blooms 


cael wits “dx =~ b, as V—->oo. 
PS Jo  \SP 


The validity of the inequalities (i.) and (ii.) depends on the form of the 
functions g(E) and h(E, V), and a general discussion is impracticable. 
Husimi (1939) has given expressions for g(E) and h(E, V) for a perfect 
gas in a cubic box, but it is more convenient, in this case, to evaluate 
the b, by the following method. If the gas is in a cubic box of side a, the 
energy levels are 


EL an, nek T (2+ m2-+n?), 


where «=h?/8mkTa? and J, m, n are positive integers. Therefore 


co 3 
d'w,; exp (—spe,) = E exp (— si) | : 
i t=1 


Now z exp (—sal?)<| exp hiner) AS nee (—sal?), 
l=1 0 = 
ee 2 Ryd ae 
and therefore 3 (=) —1 oe exp (—sal?) << 5 (=) é 
' N (2amkT\3? 
Hence 2 w exp (—sue,)< 3 | ——— for all s, 
: ph s 
and 


—dy) for N>Mands<sy, 
8 


N 
yates P (—spe;)> = 


2amkT (M\2/8 5 
where oN ae 0), 4 p OM: 


ae A 


‘Thus, we see that postulates (i) and (ii) are justified in this case with 


QnmkT\2/2 1 
irs ee pak" 
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Postulate (iii.), that 26,78 has radius of convergence unity, is not 


essential and has been made for simplicity only. If we assume, more 
generally, that the radius of convergence is 7, the condition for (19a) 
to hold becomes €)<r and (19 6) becomes 
f=—kl 2 byr?. 
s=1 

Since exp (—se;) is a decreasing function of s, we see from (5) and 
(6, i. and ii.) that 7 cannot be less than one. + will be greater than one 
only if 6, contains a factor exp (—sf), and this can only happen if 
g(E) is zero for 0<E<f/u. This case, however, can be reduced to the 
case of r=1 by a change of the zero of energy. 

Postulate (iv.) has also been made for simplicity only. If 2Xsb, is 


& 
divergent, condition (a) of (14) will always be violated before condition (0). 
This means that the gas will not condense, since equation (19b) will 
never apply. An example of this behaviour is to be found in the two- 
dimensional Bose-Einstein gas, for which b,=2amkT/ps?h?. 
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ABSTRACT. 


In an earlier paper formule were given for the macroscopic behaviour 
of a perfect Bose-Hinstein gas. In this paper formule are given for the 
mean occupation numbers of the energy levels and for the fluctuations in 
these numbers. These formule are in agreement with those obtained by 
earlier (invalid) methods. 


§1. IrTRODUCTION. 


As in the earlier paper (Fraser 1951), we shall consider a perfect Bose— 
Hinstein gas of N particles in which the energy levels of the individual 
particles are <9, €,, €2,... With weights wo, w,, w.,... We shall find. 
asymptotic expressions for n;, the mean occupation number of the ith 
level, in the limit as N tends to infinity. 

Schubert (1946) and Dingle (1949) have pointed out that the derivation 
of the mean occupation numbers by the method of steepest descents is 
invalid at low temperatures. Expressions for the occupation numbers in 
the most probable configuration of the system may be obtained by other 
methods, but these are also open to criticism (see, for example, Dingle 
1949.) It is well known that expressions for the occupation numbers 
which are obtained by the ‘‘ mean value ”’ and “‘ maximum probability ” 
methods are, in general, identical, and it is of particular interest to verify 
this equivalence of a Bose-Einstein gas, since the fluctuations in the 
occupation numbers are unusually large. The formula usually given for the 
fluctuation is 


i= (np oe are re, teste ee (TL) 


(n,;—1;)P? =n, 
and (n;—n, P=, —(%)*, 


for a Boltzmann and a Fermi—Dirac gas respectively. The validity of 


equation (1) will be discussed in § 6. 
oS Seale pt i aa lt i a ae De in a IT ae ae 
* Communicated by the Author. 
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§2. Exact ExprEssIONS FOR THE MpAN OccuPATION NUMBERS. 
We start from the partition function Zy, defined by the equation 


Lng=2) XP (— ae nei) : 
i=0 


in which p»=1/kT, and X’ means the sum over all values of the set of 


positive integers 7, 21, 2, ... such that 2n,=N. Zy may be considered 
as the coefficient of CX in : 


EZq(%— I [1—Lexp (—pe))I-*- oe ans 
N j=0 


Expressons for the mean values of the occupation numbers and their 
products can then be expressed as functions of Zy (Schroedinger 1945). 


lone 


n=—->- a etl 3 th thesis 
Nn; pLy de; ? ( ) 

Ll ®Zy 
n= —— —, ee pale ae 
as pZLy de ,0e; ( ) 


If we differentiate equation (2) with respect to «; and pick out the 

coefficient of £X, we see that (3) takes the form 
N 
n(N)=o,; 2 “N= exp (—rue)) itt. ee fs ae 
r=1 N 

Here n,(N) has been written instead of n;, since we will consider the 
behaviour of 7; as a function of the total number of particles N. 

Since Yn,(N)=N, equation (5) leads to 


N 
NZy= 2 Zy_,2w,exp(—rpe).° 2. . S «os (6) 
r=1 a 
Setting N=1, 2, 3, ... in turn in equation (6), we have, in principle, a 


method of finding Zy. Thus equations (5) and (6) are sufficient to deter- 
mine 7,(N). 
From (5) we can also obtain the important relation 


— Z 

n(N+1)= 7 (w,+n,N))exp(—pe). . . « © (7) 
N+1 

This relation has been obtained by Ansbacher and Ehrenberg (1949) who 

point out that, if we assume that n,(N-+1)=(1+-5)n,(N) and nelgect 5, we 

immediately obtain 


n; N = ed i ae a ° e r e ° e 8 
1. Calon) OD eI oe 
However, this procedure is only justified if 5 is small compared with 
1—exp (—pe,)Zy/Zy+,. In fact, if i=0 and the gas is condensed, both 
these quantities are of order 1/N and 8 cannot be neglected. For i-40 


it seems that the neglect of 5 is usually justified, but we shall make no 
assumption of this nature. 


te ee 


Condensation of a Perfect Bose-Einstein Gas.—II 167 


§3.° PROPERTIES OF THE PaRTrITiIoN FUNCTION. 
We shall define a quantity A,(N) by the relation 


Bet Nvcxhnc Dec acaral irae (mene G2) 


A,(N) will, in general, be a function of the temperature as well as of ¢ and N, 
but we shall only be concerned with its variation with N at constant 
temperature. 

Substituting (9) into equation (7), we find 


Z 


n(N+1)= 7 Az (N)n,(N) peyote: Wh RTs THO): 
N+1 
=] eal Zy41 
and Az (N+1)= +exp (—pe,) 1—)0N) eal La 
N 
Summing equation (10) over all 7, we obtain 
Z ca 
Nee A Nn iN)a te}... (12) 
Zy41 
We can now prove that 
: esicxiniiics ar forvalleN see... 1(13) 
Zy41 
Let M be the lowest value of N for which (13) is untrue, so that 
Z 
ae > exp ("Eq) eet ee ace Peg (14) 
M+1 


(equation (6) shows that M is greater than zero). The inequality (13): 
is true for all N<M and hence, by (5) and (9), A,(N)<exp(pe9) for N<M. 
Therefore, from (10) and (14),7,(M-+1)>7,(M) for all 7 and, from (7) and (14), 
M(M-+1)>1-+7,(M). Thus an; Ea ee (M), which is impossible. 


Therefore there is no value M for which (14) is tate and hence (13) is true for: 


all N. 
_ If we use the inequality (13) in the equation (5), we find 


n(N)<w,; z exp [—7u(€,—€9) ], 


rT=1 


earl dea ala . . ° ° . e (15): 


Ww 
Therefore, if Ne ee ee 2 LG) 
e540 exp [u(€; -€o)] -1 
is finite, we must have n,(N) >N—N, whenever N>N,. Thus condensation. 
must occur whenever the energy levels are such that the series defining N, 
converges, since, when the total number of particles exceeds N,, a finite 
proportion of the particles must be in the lowest state. 


€ 
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Since (16) can be considered as defining a characteristic temperature T,, 
associated with a given number of particles, we can say that a gas containing 
a fixed number of particles must be condensed at temperatures less than 
T,. It should be emphasized that this result, in common with all those 
obtained in this section, is true for a general set of energy levels (provided 
only that the series in (16) converges) and for a finite number of particles. 
In order to show that the gas is not condensed at temperatures above T, 
_and to find the exact distribution of particles among the energy levels, we 
must consider the limiting case of large numbers of particles, but finite 
. densities. 
We shall first prove some other useful inequalities which are always true. 
From equations (10) and (11) we obtain the relation 


n{N-+2)= 7 +exp (=e) (1A) J) bN+1) . (17a) 
N N 


Zy+9 


_ Zyay {Ar N+ a exp (—sL1e;) (x10) 7)} nN). 


— Lyy. N+1 D 
(17b) 
We can now prove, by induction, that 
Z 
Z AAs 1(N)S-1 forall N ands. =) a 
N+] 


We assume (18) to be true for all 7 when N=M. Therefore the inner 
bracket expressions on the right-hand sides of (17a) and (176) are positive 
for N=M._ If we also use the inequality (13) we can write 


Zyv1 Syria — — Z Z, 
po ne (MT) nM oe ee 2) ee 
Zy42 Zy are a a at An Zy } maw. 

Hence, summing over all 7 and using (12), we find 
Zyii1 Zui Z Z 
pa | a HY) fee | <M 2 M+1 ; M+1 < 
Zyr2 Ly ane st a Zyuig Ly (M2): - + (19) 


From equation (11) and the right-hand side of the inequality (19) we obtain 


V6, ig 
7 Ar 1(M+1)> Zyi4 : Zy+4 1; 
M+2 u+2 Zu 
Therefore (18) is true for N=M+1. It is true for N=0 and hence is 
always true. Consequently (19) is true for all M and can be rewritten in 
the form : 
Ly +1 N Guy Z N Zy+4 
7, N+2 a 7, : : -) omens (20) 
nee N+ nit 4y+e 5 
This relation, expressing the fact that Zy/Zy.,isa slowly increasing function 


of N, is of great importance in evaluating the asymptotic values of the 
mean occupation numbers. 


a 


7 
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§4. THE Mean Occupation NUMBERS ABOVE THE CONDENSATION 
TEMPERATURE. 

In this and the following sections we shall be conberned with the 
asymptotic behaviour of quantities such as 7,(N) when the number of 
particles N and the volume V tend to infinity in such a way that the density 
p=N/V remains finite and constant. The volume enters into the expressions 
implicitly, since the distribution of the energy levels depends on it. 


Zyey _ Zy_ Se) Ly ptt Zy_4 
Zy Zee r+] Ly pa. Ot Zy , 


the inequality (20) shows that 
Zyi\" _ Zy_, Zy_1\" 
CF) alee 
_N(N=1)....(N=r+1) 
ve Npat 


Since 


where a 
By the use of Stirling’s expansion for log (N!), we find that 


r2 r y4 
O=1— 5 + oe 4 0 (xs): 


From equation (5), we see that if M is any integer less than N, 


N Wi, : 
(—me)to, 2 [= exp (—r«;), 
=] Ly r=M+1 N 


n(N)= 


and, by (21), 


M r M 7, r 
», at exp (—7pe;)> 2 Ene Ces z a .( a =!) exp (—7pe,). 
g=1 Ly r=1 Ly Zy 


Let M=N*, where « is less than 3; then a,~1 as N tends to infinity, as 
long as r is not greater than M. Hence 


N& NOT, r 
Py ABS (—rpe)~ 2 *-1) exp (—7pe¢,). 
r= Lb =1\ 4y / 


r= NX+1 N 


NO /7, Sy ah 
(Smo, 2 (BAY exp (— nate; 2 Fetexp (me) 
N 


lee} r 7. 
—o{ a (=) exp(—me)— 2 =" exp (—me) 


r=NO+1 


foes W; meeps ; 
As ico) (EnlZa-)) xD Ged bela 
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where, by (21), ‘i 
one 
O=— RUN ce) = ee 7 ) exP (—7pe;). 


r=NO+1 N 
Thus 
Na 
a= NT yo perl tiv te lO eres lee N-1 ieee _ (22 
™u(N)~ TT a) exp (14)) at -| 52 ale (1) | } sluisteuuae 


where 0<@<1, and «@ is any positive number less than $. Since, by (13), 
Zy-1/Zy< XP (Hey), 1—(Zy_1/Zx) EXP (—p€q) is positive and less than 1. 
Thus if we write z=1—(Zy_,/Zy) exp (—p€9), 


No 
Lim Ca exp (—n«i)) = Lim {exp [—pN*(e;—€9)](1 — a) N%aley 
N Noo 


N>o 
 < Lim {exp [—pN*(e,—e9)] exp [—Nz]}. 
N-o 

Hence, if either (e;—¢,) or x is of order greater than N~, the coefficient of @ 
in (22) tends to zero as N tends to infinity. 

sa 
(Zy/Zy_1) eXp (we;) —1” 
(a) at all temperatures if («;—«))~1=o(N%), 


r (b) for all a if [1—(Zy_,/Zy) exp (—peo)|-1=0(N%), . (24) 


TNS  a(Njre (23) 


where « is any number less than $. 

(23) is, of course, the well- known formula for the occupation numbers, 
and (24a) shows that it is always true except for some levels with energy 
very close to that of the lowest level. 

We will show, however, that the condition (b) under which (23) is valid 
is satisfied for all temperatures except those in the immediate neighbour- 
hood of or below the condensation temperature. This conclusion is to 
be expected from the fact that (b) demands only that the number of particles 
in the lowest state should be of order less than N}. 

For convenience we write A=Zy_,/Zy, so that (23) takes the form 


Ww; 


A~1 exp (pe;) —1° 
So long as the condition (246) is obeyed, the identity eur (N)=N takes the 
limiting form F(A, T)=p, where 


W; 
man ik WF Texp rosie ? 


The omission of the term Wed te to 1=0 makes no difference for 
large N, since 7)(N) is of order less than N?. The equation F(A, T)=p may 
be used to define A as a function of T (for given p) and the condensation 
temperature may be defined as the value of T, T, say, for whichA=exp (pe,). 


n(N)= 


es 


Condensation of a Perfect Bose-Einstein Gas.—II 171 


This value of A is clearly outside the range of validity of (24b) and this 
definition may be considered at present as purely conventional, although 
it should be observed that it corresponds to the deanicion of the 
condensation point given, for finite N, in equation (16). Any temperature 
T’>T,, such that (T’—T,) is some finite temperature difference, however 
small, will correspond to a value of A, \’ say, which differs from the limiting 
value exp(e) by some finite amount. It is clear that, for large enough N, 
(246) will be satisfied if Zy_,/Zy is equal to )’ and hence (23) will be true at 


the temperature T’. It follows that (23) is true, in the limit of large N, for 
all temperatures greater than T.. Bos 


§5. THE Mzan Occupation NUMBERS BELOW THE CONDENSATION 
TEMPERATURE. 


In the last section we defined a temperature T, by the equation 


1 w; 
ah Ne iS exp [4 ,(€;--€o)]—1 mi, ae Boe eae 
where ,=1/kT,; we have already remarked that this is equivalent to the 
limiting form of equation (16) in which we considered the temperature as. 
fixed and the number of particles as our independent variable. It is 
now convenient to return to that view-point and to define the critical 
number of particles at a temperature T by the equation 


Ww; 
N= & ————.,,_ ve DT. 
° ¢z0 XP [M(€; <9) ]—1 Lg 
From inequality (15) we see that )(N)>N—N,, and from equation (5) 
and inequality (21) we see that 


(62) 

RIN 

"0N) < Ty] a) SP (He) 

eee vie Os aN SN! 
Oe aad Unlina) expel 
We thus obtain the inequality 

1 iors Genel tee NY (26) 
— WE P (HE N_N, © 


where the left-hand side is a restatement of (13). 
Let us consider the number of particles, N, to be such that N—N,=rN, 


where r is some positive number. Then (26) becomes 


Zy oy. 
l< Tet exp (Weg) <1 +- rN’ 


Tee ate 
similarly 1< Tae. eXP ({uep)<1+ aye 
N2 
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where M<rN. By forming the product of these inequalities with M 
running from zero to Né, where B is some number less than 1 (so that 
N’ <¢N for large N), we obtain, after some rearrangement 


exp (e,N®) > 2x—nb >(1+ read ia exp (ue,N*), (N?<rN). 
Mo TIN rN —N? o : 


If we now use this inequality for Zy_x¢/Zy in equation (5) and proceed ina 
manner similar to that used in the derivation of (22), we find 


ns sa Bileee 
nN) ~ exp lee) =I {1—6 exp [—N*u(e;— 9) ]}, 
where 0<0<1. 

If (€;—€,) is of order N-” where y<1, B can be chosen so that y<B<1 
and the coefficient of § is exponentially small for 740. For a perfect 
three-dimensional gas, for instance, y=2/3, so that this condition is met. 
‘Thus 

exp [e(e; €o)] (27) 


Ww; 


So 
iz0 &XP [w(e,—€p)]—1 
(27) has been proved for N—N,=rN, where r is any positive number, 
however small, so long as it remains finite as N tends to infinity. 
If the particle density p is such that the condensation temperature, 


defined by (25) is T,, then any temperature T’’<T, will correspond to a 
critical density p’’, where 


and hence No(N) ~N— 


1 w; 


Lt <5 & —— =", 
be fy izo xp [wu ea} p 


and p’=1/kT’’.  p’’ will be less than p by some finite difference, rp say. 


Therefore at the temperature T’’, N—N,=rN and the argument used above 
is valid. Hence equation (27) is valid at all temperatures T<T,,. 

Thus we may sum up the results of this and the preceding section as 
follows. 

If T, is defined by the equation 


1 w; 
Lt <5 & —————_ > =p= Lt NV 
ty i~0 OXP eer . V,N>0 Me 
where p,=1/kT., then 


(Zy/Zy_4) exp (ye;)—1 or 1 > te 
and 
n; N min ee F - 0 
a exp[n(e—e@)l—1? 7? 


- Cay Sd Bead be 
ny(N)=N— 2 nN), 2 


i#0 
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In obtaining these expressions we have assumed that the limit used in 
defining T\, exists and that («;—¢ ) is of order greater than N-’, yal, 


'  §6. Frucruations. 
We have seen in equation (4) that 


= OE, 
: eT. de,0e; 


By differentiating (2) twice we find 
X Zy_, exp [=(r— nei] —exp [—(r—1)ne)] 


N;N;=W,W 


exp [—H(e,+¢;)] 


gy pao. Zn exp (—pe;)—exp (—pe;) 
1 aa —— . . . 
~ exp (=e) exp (pe) + exp (—pe;)n,—w, exp (—pe)n;} if tAj, 
(28) 
and when i=), we find 
=) ae NZ, 
m= ;-o(w;+1) 2 =" (r—1) exp (—rHe), . . . (29) 


r=2 N 


(28) and (29) are exact formule. 
By arguments similar to those which lead to equations (23) and (27), we 
obtain from (29) 


HP 


(n;—n;)P?~ 


W; 

for all2 when T>T, and for 140 when T<T,. 

To find the value of n? for T<T,, we observe that 

Nnj=n2+Zingn,, 
hence - (Mp—M%)2=% X 1;— XY MN; 
; i#0 = 4 #0 
- and using (28) with 7=0, we obtain 
Ww; \ wy Oxp [1(e—<) 10; ; 
izo OXP [H(<;—€)]—1 
(30) 


=i), Ce 4 i SS 
(o ; Mo) =o 12 ms izo &Xp [m(e;—€p) ]—1 


Although 7, is given asymptotically by (27) for T<T,, we cannot necesarily 
neglect the first term in (30), since 


Ss, (7 exp [p(e—€)]—1/’ 
while neglible compared with n) X n; need not be negligible in comparison 
i40 
with the second term in (30). Since, by (5) and (18), 


Ww; 
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we can deduce from (30) that 


s Ma Pea yr nj)? exp (H<;) 
Gee aybe 5 ee eD Leer 6) Pi ye OR oo ce eee 
(Mo— No) Sori, exp [“(<,—€o)]—1 eXP (Hé9) 540 ol! 


Thus we have shown that 


—— 8 (Mi) eralla ial ok 
(n;—N,) +n, for alli if T>T,, 

w; 

roe 
(n;—n,)? ~ hea +n,;, for all 70, 

W; a 

rbd Reseed be) 
an ey: (n,)? exp (H«;) ’ 
9 8 eXP (H€9)"iz0 W; : 


The expressions for these fluctuatior.s which are obtained by the method 
of steepest descents differ from these only in that the “ less than or equals *” 


sign in the expression for (%—7)* when T<'T, is replaced by an “equals” 
sign. The expressions for i=0 are sometimes quoted as holding generally, 
even for 7=0 in the condensed state. That this cannot be so can easily be 
seen, for if 
(M%—No)?= (Mo)? + No 
(we have taken w)=1 for simplicity), then 
N}=2(Np)? +N. 

At a sufficiently low temperature more than 1/4/2 of the total number of 
particles will be in the lowest state. Therefore 2(n,)? will be greater than N? 
and hence n? will be greater than N?. This is impossible, since n»<N. 


- The formule which Fowler (1936) obtains by the method of steepest- 


descents are not open to this criticism and reduce to the formule given 
above with the inequality replaced by an equality. 


§ 7. CONCLUSIONS. 


The formulz obtained in this paper show that the expressions obtained 
by the method of steepest descents are correct, except possibly for a very 
small range of temperature near the condensation temperature which 
decreases to zero as the number of particles and the volume increase to 
infinity. The only new assumption which we have made is that the 
separation between the lowest pair of energy levels, (¢,—e9), is of order 
greater than N~’, where y<1. If this assumption were not made, 
complications would arise in the condensation phenomenon even if the 
orthodox formule were valid everywhere, since condensation would take 
place into more than one level. If, for instance, ¢,—e,=e/N, then if 


1 


"WT exp (ue)—1 


——— 


. 
OO ee eee 2 eee ee eS 
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and n)=aN, A would have the value 


Nae: (14€9) 
= 1-+-1/aN 
and n, the value 
ae N 
* pe+lfa” 


In other words, a comparable number of particles would have condensed 
into the level above the lowest. 

From the results obtained in §§4 and 5 we can deduce that the 
macroscopic behaviour of a perfect Bose—Einstein gas is identical with that 
found by London (1938) and Fowler and Jones (1938). This is in agree- 
ment with the conclusions reached in paper I. (Fraser 1951). 

The only work which appears to be in positive disagreement with the 
results obtained by the method of steepest descents is that of Temperley 
(1949). Since this work deals essentially with a finite (although large) 
number of particles, it is difficult to compare with the conclusions of 
§§4 and 5 of this paper. Temperley’s conclusion that, in his particular 
model, condensation only takes place at a temperature much lower than that 
predicted by the method of steepest descents is, however, in disagreement 
with the consequences of equations (15) and (16). This disagreement arises 
from the choice of a criterion for condensation. Temperley has required 
that, when the gas is condensed, all particles added to the gas should go 
into the lowest state. This, however, can never occur for a gas consisting 
of a finite number of particles. In fact equations (10) and (18) of this 
paper show that the addition of a single particle leads to an increase of 
all the occupation numbers. It is only in the limit of infinite volume that 
we can speak of the occupation numbers of the excited states being 
independent of the total number of particles when the gas is condensed. 
Relations deduced by Temperley for his model can be used to prove that, 
when the temperature is below the orthodox condensation temperature, 
almost all particles added to the system condense into the lowest state. 
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XVI. The Unsteady Flow of Viscous Incompressible Fluid Inside an 
Infinite Channel. 
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SUMMARY. 


The problem considered is the motion of viscous incompressible fluid. 
inside an infinite wedge-shaped channel of angle z/n, the channel being 
started suddenly from rest with uniform velocity in the direction of its. 
length. The solution is obtained by considering a suitable Green’s. 
function for the heat conduction equation. The velocity is found to be 
expressible either in terms of error functions or in the form of integrals. 
which are suitable for computation, depending on the value of n. In all 
cases the skin friction appears as a simple expression involving well known. 
functions. 

Particular reference is made to the right angled channel (n=2), the results. 
being extended, by Rayleigh’s hypothesis, to include some discussion of the 
steady flow inside a right angled channel having leading edges normal to the: 
incident stream. From this, also, an approximation is derived for the 
frictional force on the interior of a finite rectangular channel of suitable 
dimensions. 


1. Statement of Problem. Choosing cylindrical coordinates (r, 6, z), two: 
semi-infinite plates coinciding with the planes 0=0, 0=7/n, respectively, 
form an infinite wedge-shaped channel of angle z/n (n is an integer>1). 
The region inside the channel is filled with viscous incompressible fluid, 
and initially the channel and the fluid are at rest. At ¢=0 the channel is 
set in motion impulsively with uniform velocity W parallel to Oz; the 
problem is to determine the subsequent motion of the fluid. 


2. The Equation of Motion and Boundary Conditions. Let (u, v, w) be 
the components of velocity of the fluid in the given coordinate system; in 
particular, w is the component parallel to the axis =0. Then it is clear 
that there is a solution of the hydrodynamical equations, and the solution 
of the above problem, in which w=v=0, provided w(r, 6, t) can be deter- 
mined to satisfy 
Ow (dw ldw , 1 dw 
an (Ge Tor 3) . . . . . (2.1) 
a re ee ee ee SS ee ee 


t Communicated by Dr. A. Fletcher. 
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and the boundary conditions 
w=0 when ft=0, 
w=W if 6=0, 6=n/n, 

dio Vig alae sre gain | When.tm0, <j). : (2.2) 

(v is the coefficient of kinematic viscosity.) 

3. The Solution of (2.1) It is convenient to write 

Te | ES ren re et. ae (8.1): 


so that w* must also satisfy the equation (2.1) and the boundary conditions. 


Ce Vs when 7_-0; 


w*=0 if 6=0, b-=n/n | 
when ¢>0. - . (3,2) 
w*——W as r—> (940, a/n), | 


The expression for w* can be obtained directly by considering a Green’s. 
function relating to the equation (2.1) and appropriate to the boundary 
considered (Carslaw and Jaeger 1947). An expression is obtained for 
w*(r, 0, t) at the point (r’, 0’), say, and the required Green’s function, wu, is 
of the form 

eee NOt) Gig t=) ew ees. te. Dee (As) 


(in this, and in all that follows, we take t>7.) 
We require u to satisfy the equation 
Ou (aw low 1 ow 
a \Or | ror | 7 Oe)’ 


and the following conditions 


Gi) 


ID Saye reed OE Mia i alt i Aik veel ee OLEES 
t+>7r+0 


everywhere in the region S(0<0<z/n), except at the point (7’, 6’) where 
' 1 1 2 19 y y s 
+ — —___ r’2— 2rr’ cos(d@—O » | 5 gw (ivy) 
4 4arv(t— 7) exp| 4v(t—7) an ( 
and also u=0 if 6=0, A=n/n. Seu et see ee ae 
In addition, we shall require u and @u/06 to tend to zero in a suitable manner 
as 7——> oo , to ensure the vanishing of a certain line integral. The function. 


u which we shall choose will certainly satisfy this last condition. Finally, 
under these conditions, and the conditions imposed on w*, the required. 


expression for w* is 
SAL 00 4 ; 
17 Ue )=W | | Gf adds «ae a 2 Kes 
‘ 0=0/ r=0 


It remains, therefore to select the appropriate function w. 
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Consider 


1 
Ur, 9, 7, 8’, t—7T)= 


4rrv(t—r) 


GB exp— ms fae —2rr’ cos (0 _ {e+ a \ 
a 
AC, —exp— a(t moat —2rr’ cos (0 — {2¢m-41)2 1l)-—— +). 


(3.4) 


‘Then wu satisfies (ii.), since each term is a solution of (ii.), and clearly all 
the remaining conditions (iii.)-(v.) are also satisfied. Moreover, w and 
-0u/00 —-+0 suitably as r—— oo, so that 


oe ioe) 
w*(r’, 0, )=W | ; | u(r, 0, 1, 0", tyr dr db, 
6=0/ r=0 


x 


or, reverting to the usual notation, 


mln ee) 
w*i(r, 0, )=W | | a2; a, 7;:0;. te dziday i. Fr. a0) 
a=0/4 z=0 


where 
1 


4nvt 


i Pee pa i 2m 
or exp— bi)” +7?— 2rz cos arr ee 
ne oat 


om ot eae De 
ee a be 


It is at once evident that the expression (3.5) satisfies the boundary 
conditions on 6=0, 6=7/n, since the integrand w is zero for these values of 0. 
The expression can be reduced further by carrying out the integration with 
respect to z, and, indeed, except for a certain group of the integers n, the 
integration with respect to « can also be performed, the expression then 
reducing to a (finite) sum of tabulated functions. 

Writing p?=1/4t, c=cos (@—B), s=sin (@—f), we have 


al _.° ; xp{— gat 2re00s{0—B})} ads 


u(z, %, 7,8, t)= 


=p" exp (— prs?) | exp {—p?(z—rc)* }zdz 
z=0 \ 


=p” exp (—p?r?s) | exp (— px) . (atre) dx 


= fC 


+ 
=} exp (—°)-+y¢. exp (—7?s*) 5 {1+ erf (0)}, 


“where n= ’ 
7 


r Heer ee 2 i x2) d 
Zp? er W)=5 | exp(—2") op 


~ sre 
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Hence, placing y,,=0—2mz/n 


W aln n=1 
pees 
acl, Bs 
| © €08 (Ym_—04) . EXP (—7? sin? {y,,—« }){1 + erf (7 C08 {/,—c})} . 
cS Oo 
— 608 (7m 4.1 +a) .exp (— 7? Sin? {Ym 41+ o}) {1+ ert (4 COS {741+ «})} 
(3.7) 
or, writing f(y, b)=cos ¢ .exp (—7? sin? ¢) {1+ erf(ncos¢)}, . . (3.8) 
ua ~ (2m-+1)= —6 6—(2m+1)= 
wma 2 Llane, Worsae [oo tania syas 
_ m+ 1)— 
ais ie 1 p@m+1)2 —6 2(m-+1) = —@ 
cal Fad) ab— fn" aftr d)48 
eae 0 2m = —6 (2m+1)> —8 
(since f(7, ¢) is an even function of ¢) 
W 2n-1 . (m+1) —6 i 
<5 2 (Um | , ” alr dds. 
we ] 20-1 E r(m+1)= ~6 F Ar 
Thus Wa 2, (—) lee. nf(ad)dd —« .* . (3.9) 
| 2n-1 /n mir 
tae Da —])m ——@)d¢d. . (3.10 
at 2 (—U" [at (n b+ 52-8) ap, . (8.10) 


From these last results it can easily be verified that w* does, in fact, 
satisfy the remaining boundary conditions (3.2), the first and third of these 
now being included in the single condition w*——> W as y—> «0 (0<@<m/n). 
For convenience, take n>2, and consider 


where bm =mr/n—O. 


Then, since 0<0<7/n, the range ¢,,<$<bm+7/” excludes the points 
$=0, d=7 except when m=0, n respectively. Thus, if mA0, » 


a es *ncosg.exp(—n2sin’d) .{1-+ erf (700s $)}d¢ 
= exp (—7? sin? €)7 cose [ 3 {1+ erf (7 cos ¢) }d¢, 


where ¢,,<&<$,1+7/n. 
Hence, since siné +0, and Jerfa]<1, 


Lj—0 48) 7 0. 
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In the cases m=0 and n, cos ¢ does not change sign in the range of 
integration, so that 


I,=(1+ erf (neos€)} [" 


=" {1+ erf (7 cos &)} ext (n sin (4m =) ) — erf (7 sin bn) | , 


Now, if m=n, cosé<0, so that erf (7 cos £) >— 1 as 7 > 00, and therefore: 
I,—>0. But if m=0, cosé>0, and ¢).<0<¢,)+7/n, and therefore: 
I,—~2n! asn—+oo. This establishes the required result. 

Thus, finally, the expression for w is given by 


mu 
n 


ncos¢.exp (—7? sin? ¢) dd 


nf(n,¢)db,- . » (3.11) 


w ] 2-1 2 (m+1)= —6 
Tm ea) 
Henceforward we shall consider the simplification of this result when n 
takes certain values. Accordingly, it is convenient to introduce the 
siffix n, w now being replaced by w,, and so on for other quantities we shall 
consider. In particular w, and 7, will refer to the velocity and skin. 
friction respectively for the standard flat plate problem of Rayleigh (1911). 


m-—6 
n 


Writing g(y, $)= cosd.exp(—y*sin?¢d) . . . . . (3.12) 8 
and h(n, 6)=cos .exp ae sin? ¢)erf(ncos¢), . . . (3.13). 
then Fn, P)=9(n, $) +h, $); 
and pasta =—g(n-. ¢), : 


A(n, 6+7)=h(y, $). 
The expression (3.11) now simplifies as follows :— 
(i.) n=2p-+1 (p an integer>0). The terms involving h in (3.11) vanish,. 
and hence 


Wp 1 2 (m+1)= —6 
Weed oy mike x n cos fd. exp (—7? sin? J) dd 
a ‘ 2p ; 7 . [mr 
=1—} 2 (—1)™| erf( ysin< (m+1)-— — —erf (sin {= ob) |, 
m=0 n n 
i.e. m14 FS (—1)ert( naind™ 9 3.14 
<i ee. n mA signs? wees Se 


(ii.) n=2p (pan integer>1). In this case the terms in g in (3.11) vanish, 
and we have 


Wn 1 "PS (m+1)= ~6 
Pe Pi : 
W a ee I le % 7 cos ¢ .exp (—y* sin? 4) . erf (71 cos 4) dd, 


(3.15): 
(ii.a) There is a further simplification of 0. 15) if n=4p+-2(p an integer>0).. 
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For, writing 


ba (m+1)= —6 
In=(—1) per 7 COS $ . exp (—7? sin? ¢) . erf (7 cos 4) dd, 


then 


+(m+1)= 6 
AapsS 7 cos ¢ .exp (—7* sin? ¢) . erf (4 cos ¢) dd 


hel 


2p+1+m— 


Ed 
(m+1) A —6 


e=(-—1 jet? ee sin ¢ .exp (—7? cos? ¢) . erf (1 sind) dd, 


and therefore 
Dt Tep stam 


ee 
ao Sh 


lear 1)= -6 n COS p . exp (—7? sin? 4) . erf (7 cos $) 
: d 
—7sin ¢ . exp (—7? cos? ¢) . erf (y sin d) 
(m+1)= —6 


Tt 
oa a i a [ext (1 cos ¢) . erf (7 sin ¢) iE 
Thus, finally 


2 
7 =1+ 2 (ment (1 003 | —o ‘ ext (y sin -0 ) a WED) 


The results (3.14) and (3.16) give simple expressions for w, in terms of 
error functions, and deal with all cases of n with the exception of the group 
n=4p (paninteger). In this case w, may be computed from (3.15). The 
first three examples of w,, are 


w,=W[1—erf(nsin@)] (Rayleigh’s case), 


—6 


W,= W[1—erf (7 cos @) . erf (7 sin 4)], 


e 1—erf (7 sin @)—erf (7 sin {7/3—6}) 
erf(ysin {27/3—6}). 

In all cases the velocity distribution must be symmetrical about the plane 
6=0,—=n/2n. Fig. 1 shows 1—w,/W plotted against #,—@ for n=1, 2, 3 
respectively in the cases n=2 and n=3. 

4, The Skin Friction. Since the motion is symmetrical about 6=6,, it 
is clearly sufficient to consider the skin friction 7, on one plate only. 
Selecting the plate 0=0, then for t>0 


(3.17) 


Ws; = 


ee OW, i KGLP = 
(a Sp met Ek, (evaluated for 6=0) 


a “E (-1) ais abate mm) s(n ere *) | 


i Wikws (ye thay ae sage ee LT 
or en ee (x, =). (4.1 
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(4.1) gives an expression for 7,, for all values of n, and it follows immediately 
that 

Ty —>T, 88 so. «wwe ws (42) 
That is, sufficiently far from the corner, the skin friction approximates to. 
Rayleigh’s value for the flat plate problem. 


For, if m<0,n, then sin m/n40, and f(y, mn/n)>0 as n> 


while if m=n, f(n, m7/n)=—(1—erf (n))—0 as 7-00 
and if m=0 f(n, ma/n)=1+erf (n)>2 as ->0O. 
If nm>1, then 7,=0 when re re ee TEE 

2n—1 n=1 

since Te (O) == aad py (—1)” cos = 4 (av an ‘i 

2(7vt)* moo n 
0 (n>1) 
Fig. 1. 


Further, if n>1, it follows that 
T1—T,>0. & oN 9 eee 
To establish this, consider the following simplified expressions for 7,— 
(i.) n=2p-+1 (p>1) 


2uW P 
T1—-T,= at (oat ee Bs 1 cos “7 -exp (=n *sine 2) . (4.5) 


ee ee 
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(il.) n=2p 
uW [ oi es 
alee ert ne > (— 1 ett — 
(arvt)? ca Paes an 
nam] xexp(—ytain™) et (o00™)] p> 
pW 
vit Lert (p=1). 


(4.6). 


Now 1—erf(n)>0, and the terms in the series (4.5) and (4.6) are 
alternately positive and negative and decreasing, each series beginning 
with a positive term. This establishes (4.4). 

Since 7;—7,—>0 as 7 ©, it is possible to determine a, such that 
7, is within 1 per cent (say) of 7, when y>«,,. 

Using (4.5) and (4.6) the following approximate values are easily 
obtained 


A a— 1 “2, 
&g—=2°5, 
oT ed oe, 
%—=3:8. 


As one would expect, «, increases with n, that is, with the sharpness. 
of the corner. 

Thus in all the cases n>1, the skin friction vanishes at the corner 
itself, but increases away from the corner until it attains the value 7, ; 
this occurs (to with 1 per cent) when 7 exceeds «,. 

Frictional Force parallel to Oz. The frictional force F,, per unit length 


6 


on the section 0<r<b is | T, ar, 80 that 
0 


Ff,—f,= ib (T1—Tn) dr=2(vt)* ie (71 —Tp) An, meee (47 


0 


, b 
where b’= 30) . 


Assuming for the moment the convergence 0 
(4.7) may be written 


f the infinite integrals, 


F,—F,,=E,—2(vt)t | (REA a sae Ce) 
2 


where . B,=2(v1)t | (z,—-T, WON. Geeta dae Woe 
0 


(so that H,>0 if m>1). 
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The convergence of the infinite integrals follows from (4.5) and (4.6) since, 
in the first place, 


5d 1 
J (1—erf n) dn=x(1—erf x)—x,(1—erf x,)+ 7 [exp (—2?)—exp (—2)], 


and as x—> oo, x(1—erf x) > 0, so that 


ie (1—erf 9) d= exp (—x u2)—a,(l—erfx,), . . (4.10) 


while the integrals arising from the terms of the series converge, since 
sin mn/n<0. 

Writing L,=yW2,, Fen+1 and A, can be found directly, and in certain 
special cases so can F,, and dA,,. From (4.1) and (4.9) it follows that 
A,, is a function of only. 

The terms on the right-hand side of (4.8), and subsequent approximations 
to them, can be calculated by considering Q,,(X), where 


__ A(vi)t (” 
Q,(X)= Ww (7 Ta) dn, 


(o) 2n— 
4. €. Qn(X)=5 | | 2+ z (—1 mer (on ~) ] dy. Soa (4.11) 
Thus A=, (00:54) 2. ne ee 


@),(X) decreases as X increases (since 7t,—7,>0), and it is clear from 
(4.5), (4.6) and (4.10) that Q(X) 0 as X > o. 
From (4.5), for n=2p+-1, we have 


Qu(X)=2 F (—1)m+1 cot (2 = =] [ 1—ext (x sin =) | oe eee 


and therefore 


A, —? Day —_}\m+1 ner 
x (—1) cot (™). 0 Ae ren A ee 


@,,(X) can also be evaluated in the special cases n=2 and n=4. Using 
(4.6) and (4.10) 


oe wes Ge 
Qu(X)= 5 | rex (—¥)—X(1—erf 2) | ; we eS 


and 
OG (Aj= = E exp (—X*)—X(l—erf X) ] 


= exp (— 47). ot ( 1) a 


2 al 
=3|s exp (— X#)—X(1-erf X) | +1—err* (x), - (4,16) 


from which A,=2/, Ay=(2/7)+1. 


— 
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A first few of the values of r,, are 


Ap=0-6366, 
Ag=1-1547, 
\,=1-6366, 
A;=2-0530, 
A, =3-2412. 


An approximation to the term on the right hand side of (4.8) can always 
be obtained. This term will be within 1 per cent of Z,, provided b’>B 
where a 

fps eg aU OA aly waa acer Py, flee eed SS (4.17) 
B,, can always be so chosen, since Q,,(X) is positive, continuous, and —> 0 
as X —-> «©. We obtain the following approximate values 


Bo=1°6, 
B3=2:1, 
By=2°6, 
B= 3-2. 
Hence, provided 6’>8,,, a convenient approximation to (4.8) is 
P= PEW, =nW | —A, | eee oa, ean (41.8) 


5. The Related Steady Flow Problem. Howarth, in his discussion of 
the corresponding problem for a semi-infinite plate (Howarth 1950) has 
used a hypothesis by Rayleigh (1911) to obtain information about the steady 
flow past an infinite quarter plane, and the same application can be made 
here. Thatis, in the results already obtained, we replace t by z/W (so that 7 
is replaced by $r(W/vz)?) and, writing w’ in place of W—w, we expect 
to derive a qualitative picture of the steady flow w’ of a stream inside the 
corner formed by the quarter planes 0=0, 0=7/n, z>0, the undisturbed 
velocity of the incident stream being (0,0, W). Of some special interest 
is the case n=2, which will be considered now in further detail. 

The Right Angled Corner (n=2). From (3.17) we have 
Perf (sett (Y) ee (Aes OGL) 
W 
writing 7 cos 0=X, 7 sind=Y. 
It follows then, that if, say, X >2 


a7 alert (¥), Wad. {th 2 AER 


approximately, since 1>erf X>0-995; that is, the flow is very nearly 
the same as that deduced by Rayleigh for the infinite flat plate. This 
implies that the disturbing effect of the corner is confined, for all practical 
purposes, to the region 0<X <2, 05 Y <2. 


SER. 7, VOL. 42, NO. 325.—FEB. I95I oO 


186 L. Sowerby on the Unsteady Flow of Viscous 


Rayleigh’s hypothesis would give, for the steady flow problem, a 
velocity pattern of the form 


U 


7: =erf (X") serf (Y’), |. - "Seeman 
W\t Fae ALA: 
where X’=he (=) ; Y’=hy (=) k 


From the point of view of boundary layer theory, therefore, we should 
expect some thickening of ‘the boundary layer in the vicinity of the 
corner, as indicated in fig. 2, while away from the corner the divergence 
from the flat plate flow will be negligible. Now the exact solution of the 
boundary layer equations for the steady flow past a flat plate (Blasius 1908), 


Fig. 2. 


Ne 


x=3 avast x’ 
Boundary layer in right angled corner. 


shows that the mainstream velocity is attained when Y’=3, and not 
Y’=2 as indicated by Rayleigh’s hypothesis. Accordingly, in this 
problem, it is probable that the influence of the corner on the boundary 
layer is effectively confined to the region 0<X'<3, 0<Y’<3. Failing an 
analytical solution of this problem, a numerical approach would require an 
assumption such as this. Carrier, in his solution (Carrier 1946) makes this 
assumption, though Howarth has pointed out that certain of the equations. 
of motion have been disregarded, and the solution is therefore open to: 
objection. 

Referring back to the unsteady problem, we have for the frictional 
force F per unit length of a section of width 6 the approximate result 


(using 4.18)) 
b 2 
oH resting | ext wc ea 
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provided 6/2(vt)t>1-6. From this we can deduce an approximation to 
_ the frictional force when a rectangular channel of sidés a and 6 is set in 
motion in the same way as the original boundary. Bearing in mind the 
requirements of the velocity distribution (5.2) then, provided 


a, b>8(vt)}, Meee Roel. (5.5) 
the frictional force per unit length on the channel would be given by 
= a+b 8 any 


taking into account all four interior faces of the channel. 
If we now extend this result, by Rayleigh’s hypothesis, to the problem 
of steady flow through the interior of a rectangular channel of finite 


length /, then if 
a,b>8 aN 
ome W jae 


the force per unit length at distance z would be 


mW | a+b) (=) —£], oo auvk SiBealtou05-8) 


and the skin friction coefficient C; for the interior of the channel would. 
be given by 


1 : W\t 8 
C= Se | [0+ (ee) 2] 


v \i 4] v \t 
=4 (=m) 1-5 (=n) it . 55 5 5 < 5 (5.9) 


The result (5.9) would be of greater interest if corresponding results 
were available for the flow on the outside of the channel. That is, it 
would be useful to extend the investigation to wedge-shaped channels of 
any angle 8. For this, a solution of (2.1) as a Fourier series in 0 would be 
available, but it should also prove fruitful to investigate a solution of period 
28 in @ given by Carslaw (1910) for use as a suitable Green’s function. 
In this latter case it may not: be possible to reduce the expression for w to 
single integrals, but nevertheless the method may yield an expression for 
7 in terms of an integral, and this would be more useful than the infinite 
series in any discussion of the limiting value of 7 remote from the corner. 


T should like to thank Professor Howarth and Mr. K. Stewartson of Bristol 
University for their kind assistance to me in the preparation of this paper. 
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XVII. Emission of Energetic Helium and Lithiwm Fragments in Nuclear 
Explosions. 


By S. O. C. SORENSEN™. 
The H.H. Wills Physical Laboratory, University of Bristol +. 


[Plates VII. & VIII.] 


SUMMARY. 

A phenomenological description is given of the processes which lead to 
the emission of energetic heavy particles from cosmic ray stars. The 
frequency, angular distributions and energy spectra of the particles have 
been determined. These, and an analysis of the associated stars, indicate 
that the particles cannot be generated by an evaporation process. On 
the other hand, direct knocking on of the heavy particle by the primary 
can be excluded. The fragments must therefore result from’ interactions 
which are intermediate between these two extreme possibilities. No 
connection has been found between the ejection of energetic heavy 
particles and meson showers. 


INTRODUCTION. 


As a result of experiments with photographic plates exposed to cosmic 
‘radiation, several reports have been given within the last few years of 
observations on energetic heavy nuclear “‘ fragments” ejected during 
the explosive disintegration of nuclei. A striking feature of these 
phenomena is that the heavy particles are frequently emitted with energies 
much greater than that which can be accounted for in terms of the electro- 
static repulsion of the residual nucleus. The problem therefore arises as 
to the physical mechanism through which it is possible to endow these 
particles with kinetic energies greatly in excess of their binding energy 
without disrupting them in the process. The emission of the particles is 
rare, and it appears likely that an extensive investigation will be necessary 
to explain their production in detail. It therefore appears appropriate 
to give a preliminary report on the material already available. As the 
expulsion of heavy charged particles presents similar problems to that 
of the ejection of «-particles, which occur much more frequently, a closer 
study of the latter has also been included in this investigation. 

The plates used in the investigation were coated with Ilford G5 emulsion, 
400 » thick, and were exposed by means of free balloons at geomagnetic 
latitude 54° N. In order to have well-defined conditions for the experiment 
only material exposed at approximately the same altitude has been 
employed. This was obtained in two level flights of which details are 
given in Table I. 


* Communicated by Professor C. F. Powell, F.R.S. 
+ On leave from the Physical Institute, University of Oslo, 
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TABLE I. 
Plate series ‘Attitude Mass of overlying Duration. of 
air exposure 
HA, 21 68,000 feet (52-+1 gm./cem.?) 100 minutes 
HA. 30 64,000 feet (60+1 gm./cm.?) 5-5 hours 


The present paper describes observations on (a) the relative numbers 
of ejected lithium and helium nuclei of great energy ; and (b) the distri- 
bution in energy of the helium particles and their directions of emission 
with respect to the “ primary ” particle producing the disintegration. 


CHARGE DETERMINATION AND RELATIVE FREQUENCY OF FRAGMENTS. 

It has long been recognized that some of the heavy particles which 
are emitted in cosmic-ray “ stars’”’ are more massive than «-particles. 
Apart, however, from the well-known lithium fragment (Franzinetti and 
Payne 1948) leading to the characteristic “hammer tracks”, it is in 
general difficult to obtain an accurate determination of their charge. In 
a recent report Perkins (1950 a) has made an extensive analysis of the 
emission of “low energy fragments’ and has discussed the various 
methods of charge determination. 

For the more energetic heavy fragments the most characteristic feature 
of the tracks is the large number of “ knock on ”’ electrons or 5-rays which 
the particles produce in their passage through the emulsion, and a study 
of the distribution of these 5-rays along the tracks provides a precise 
method for determining the charge of the fragments (Sérensen 1949). 

The present investigation has been confined to particles which are 
emitted from “stars”, and which have a residual range greater than 
700 . in the emulsion. The corresponding minimum values of the energy 
are 45 MeV. for an «-particle and 90 MeV. for Lif. Hf particles of shorter 
range are accepted, it becomes increasingly difficult to distinguish between 
these heavier particles and singly charged particles on the basis of d-ray 
counts. On the other hand, an increase of the lower limit of the range 
results in a large decrease in the total number of tracks available, and thus 
leads to results of lower statistical weight. The results of d-ray counts 
on 372 tracks is shown in fig. 1. This figure shows the relative frequency 
of occurrence of tracks with different numbers, N, of 3-rays of 3 or more 
grains, in the last 700 1 of the track. We have here included some tracks 
which have been found in other exposures, and a few tracks which entered 
and came to rest in the emulsion. ; " 

Two distinct groups of particles are immediately apparent in fig. 1, 
with a mean total number of 5-rays ~ 12 and ~ 30. Similar counts on. 
protons give a mean value of ~ 3 in the last 700 p of the track. ct 
indicates that the most pronounced group (L.) is due to belium a oh 
and the second group (II.) mostly to lithium nuclei. in the flights re ve 
to in Table I., 271 particles of charge 2e and 10 of charge 3¢ have been 
found so far, making the ratio. 

Lithium 


pO 0S eee ee eee ee ED 
Helium 
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It is reasonable to assume that the given ratio represents an upper limit 
only, since some helium tracks of unusually low 6-ray density may have 
been missed in scanning the plates. This ‘loss’? may account for the 
slight asymmetry in the distribution I. j 

That the great majority of the particles are stable or long-lived isotopes 
of the corresponding elements is proved by the fact that no secondary 
particles, such as electrons, are observed to be emitted at the point where 
the particles come to rest. None of the lithium fragments showed the 
characteristic ‘“ hammer ” tracks due to Li’. 


Fig. 1. 


Number of events. 


0 10 2 30 40 


Number of 5-rays. 
Histogram showing the distribution of the total number of 8-rays (> 3 grains) 
in the last 700 of the track of the fragment. 


In neither of the two flights have any long tracks of heavier particles 
been found although a great number of heavier ‘“‘ low energy fragments ” 
and “ taper tracks’ have been observed with a range less than 700 [Le 


ANGULAR DISTRIBUTIONS. 
In those cases in which the heavy particle has been emitted from a 
nuclear explosion of which the incident particle can be distinguished, 
the angle, y, between its track and the line of motion of the “ primary ” 
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particle has been determined. It has been assumed as a sufficient 
criterion that a track is that of a particle producing a particular nuclear 
explosion if its specific ionization is near the minimum value for charge 
Je]; and if it is the only such track in the “ upper hemisphere ” of the 
star. Camerini et al. (1949) have shown that the great majority of tracks 
so defined are due to particles which approached the nucleus and caused 
its disintegration. 

Histograms showing the distribution of the observed values of y, for 
helium particles and heavier fragments respectively, are shown in fig. 2. 
The curve gives the angular distribution of the particles with charge 2e 
in number per unit solid angle. In making these observations the spatial 


Fig. 2. 
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Value of y in degrees. 


i istributi tween the direction of 
Histogram showing the distribution of the angle be | 
on of ie primary and the heavy particle. The full cure shone tas 
results transformed to correspond to the intensity per unit solid angle. 


orientation of each track, taking into account the angle of ae in the 
emulsion has been determined. It will be seen that the ogee in 
direction is usually less than 90°, but that much larger values are sometimes 
; ae number of heavier fragments is very small, but it eae ne 
there is a collimation in the direction of the primary, a a te Ae 
with the observations made by Perkins on energetic heavy es 8. ‘ 
With the criterion for the definition of a primary peas € ee 
above, ‘‘ stars’ produced by particles ae charge ae = begs ee 
eliminated. A large fraction of the “ stars ” from which heavy fragme 
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are emitted (~ 20 per cent) are produced however, by fast «-particles at 
the altitude of the exposure and a special investigation of the angular 
distributions in these cases would be of interest. The difficulty in 
identifying the «-particle primaries except in those favourable cases in 
which the tracks are longer than 3 mm., has made the available material 
too small to allow any significant conclusions to be drawn. 


Tue DISTRIBUTION IN ENERGY OF THE HELIUM PARTICLES. 


A number of recent papers have dealt with the energy spectrum of the 
various charged particles emitted in the disintegration of nuclei with very 
high excitation energy (Harding et al. 1949, Perkins 1950 b, Page 1950). 
The main purpose of these investigations has been to check the “ evaporation 
theory ” of highly excited nuclei (Weisskopf 1937, Bagge 1944, Le Couteur 
1950) and satisfactory agreement has been obtained between theory and 
experiment. The particles considered in these investigations have energies 
up to ~ 50 MeV., but all the authors report the existence of a high energy 
tail’ due to processes which are probably totally different from that. 
of evaporation. The new method of combining “ multiple scattering ” 
and “ grain density ”’ (Fowler 1950) has made it possible to investigate 
this high energy “tail”, and energy spectra of the ejected mesons, 
protons, and deuterons have been determined (Camerini et al. 1950). It 
is reasonable to assume that most of the particles considered in the 
present investigation are the corresponding particles with charge =2e, 
and it is of interest to compare the spectra for the different groups of 
particles. 

Owing to the high grain-density in the tracks of non-relativistic particles, 


with charge =2e in the new electronsensitive plates, the energy cannot be ~ 


measured by grain counting. One way of overcoming this difficulty is 
to measure the ranges of the particles and to extrapolate range-energy 
relations, a method which is necessarily confined to particles arrested in 
the emulsion. However, an ambiguity exists in the application of 
range—energy curves, for the particles of charge 2e and 3e because of the 
presence of different isotopes. It has been shown, for example, by 
applying scattering methods to the longest tracks in the group, that 
both isotopes He*, He* contribute to the main group I. of fig. 1. It has. 
not yet been possible, however, to determine the ratio of the two types. 

Assuming all the particles in fig. 1 with less than 25 8-rays in the last 
700 w of the track to be «-particles, and applying geometrical corrections 
for loss, we obtain the differential spectrum, fig. 3, which shows the 
kinetic energy, KE, of the helium nuclei at their points of creation. 
Isotropic angular distribution with respect to the vertical was assumed. 

The longest observed track of a particle with charge 2e, ending in 
emulsion, has a range ~ 1 cm., a value which corresponds to an «-particle 
with an energy of about ~200MeV. No helium particles with 
relativistic energies have been observed ejected from stars as secondaries, 
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except in the nuclear collisions produced by heavy primaries in which they 
occur as fragments in the “ break-up ” of heavier nuclei. (Bradt and 
Peters 1950, Dainton and Kent 1950). 

A comparison of the energy spectrum in fig. 3 with the corresponding: 
results for protons, and deuterons and tritons as determined by the- 
scattering method, suggests that they are rather similar in form. The. 
spectrum may be approximately represented by a power law of the form 


const 
N(E) dE = aE 


LAD Seles Mica howe os 


where K is of the order of 1-6. 


Intensity in arbitrary units. 


40 60 80 100 [50 


Energy in MeV. 


Differential distribution in the kinetic energy of «-particles ejected from stars.. 


The lithium fragments are too few to deduce the form of their Snerey 
spectrum, but all the fragments which have been found bg end in the 
emulsion (hake a range less than 1600 p. This value corresponds to an 
energy of 145 MeV. for Lif. 
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CHARACTERISTICS OF THE “‘ STARS’ ASSOCIATED WITH THE EMISSION 
or Heavy “ FRAGMENTS ”’. 


As a second mode of attack on the problem, a statistical analysis of the 
-disintegrations from which heavy particles are emitted has been made, 
and a determination of their frequency of occurrence. Out of a total 
number of 15,000 stars in the HA. 21 plates, 189 stars of 3 or more prongs 
‘have been found from which emerge particles with charge >2e, and of range 
+700. The total number of stars giving rise to the events described 
in the present paper is 400. In this analysis we have included stars 
leading to ejection of heavy particles, which although leaving the emulsion 
before being brought to rest, have such a track-length and d-ray density 
as to indicate that their total range in the emulsion would be greater 
‘than 700. Stars caused by heavy primaries have been excluded, 
since the mechanism of disintegration in these cases may be different. 
‘Stars of the latter type commonly show a marked asymmetry due to 
collimation of the secondary fragments with regard to the line of motion 
-of the primary. Pl. VII. shows a photomicrograph of an example of this 
type, in which a heavy primary with charge Z=12+le collides with a 
nucleus of silver or bromine in the emulsion and produces an energetic 
‘disintegration with 38 heavily ionizing secondary particles. Several 
-examples of this type of event have been found in the two flights. The 
most striking feature of the event shown in PI. VII. is the expulsion of 
-a very heavy non-relativistic fragment of range 550 u. A comparison of 
this track with others due to particles which have been identified indicates 
‘that its charge was Z=5-+le. 


(a) The distribution in N,. 

Using the terminology introduced by Brown et al. (1949), we may 
-divide the tracks of particles emitted from stars into three groups 
according to their grain density g, as compared with g,,,, the grain density 
-of a relativistic particle of charge | e|: 


(1) “* Thin ” tracks gin =9<1-4G;yin (Shower particles). 

(2) Grey ” tracks 1-49 nin <9 <59min- 

(3) “ Black ” tracks g>5g,.in- 
‘The sum of the “ grey ” and ‘ black” tracks from a star is represented 
by the symbol N,, and the number of shower particles by n,. 

Fig. 4 gives the distribution in N,, for all the stars which emit energetic 


-helium particles: the observations have been confined to stars with 


n,=1 or with a singly charged relativistic primary. The corresponding 
values of the total release of energy are indicated in the figure. 

The distribution is seen to extend up to N,, ~ 35, with a broad maximum 
for N. 1117, 

The probability, P, that a star with a given N, has an associated 
‘energetic helium particle is given by a comparison of the distribution 
-shown in fig. 4 with the curve showing the frequency of occurrence of all 
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stars of this class. This is known for the HA. 21 plates (Table I.), and 
fig. 5 gives the result plotted in logarithmic scale. The figure shows that 
there is an extremely rapid rise in the probability of emission of energetic 
heavy particles with increase in N n Uf, for a given N,, the number of 
stars which emit energetic helium particles is represented by Ny,, and the 
total number of stars by N,,,,), the results of fig. 5 may be represented 
approximately by the relation 


N ON 
P—— Fe — const. e2Na, wh Oe ue So Te 
Niotat 


Fig. 4. 


Average energy release in MeV. 


300 1000 2000 


Number of stars. 


Number of heavily ionizing particles N,,. 


Histogram showing the distribution in N, of stars emitting a single energetic 
helium particle. 


The size distributions of stars emitting energetic helium particles or 
i clei is given in fig. 6. ent 
Beeieeers aitirn is distribution in N, of all the stars which et 
ane or more fast heavy particles (R > 700 < ca Re eae 
i ithi fragment. The displacement 
or one «-particle and a lithium smi ee 
whole nisin batidi towards higher N, is immediately poe . A He 
of this type has been observed in the lighter elements (C, N, OQ) in @ 
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emulsion. Pl. VIII. shows a photomicrograph of a star in which a helium 
particle and a lithium fragment are ejected, simultaneously with range 
respectively 1200 and 1000 pu. 

In order to see if the distribution in N , varies appreciably when particles 
of different energy are considered, the observations on particles with 
charge 2e were divided into two groups according to whether the range 


Intensity in arbitrary units. 


10 20 30 


Number of heavily ionizing particles N,. 


The probability of the ejection of an energetic helium particle as a 
function of N,. 


of the particle was in the interval from 700 p to 1500p or greater than 
1500 4. The corresponding values of the kinetic energy for «-particles 
are 45 MeV.<E<70 MeV. and E>70 MeV. respectively. The resulting 
distributions in N, for the associated stars are shown in fig. 6a and ¢ 
It is found that, within the statistical fluctuations, the two distri hating 
are similar in form. It appears therefore that no very strong correlation 
exists between the particular energy with which a helium particle is 
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emitted and the energy release in the associated star. In this connection 

it may be remarked that the longest «-particle observed (~ 1 em.) che 

gy ved from a star with N, only 4. Further, a number of heavy particles 
oi been en aa energies up to ~ 140 MeV. which are not associated 

ius) © other alt lair ; 4 . : 4 

ea y eavily lonizing particles—a point which will be discussed 


Fig. 6. 


Number of stars. 


Number of heavily ionizing particles.N,,. 


Histograms showing the distribution in N;, of stars emitting :— 
(a) A single helium particle with range between 700 and 1500p. 


(b) Simultaneously two heavy particles. 
(c) A single helium particle with range greater than 1500,. 
) 


(d) Energetic fragments with charge Z = 3. 


The distribution in N,, of stars which emit energetic heavy fragments 
with charge >2e is shown in fig. 6d. We have here included a few stars 
which emit particles heavier than lithium fragments and which are found 


in other exposures. 
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(b) Heavy particles and shower production. if 2 
Heisenberg (1949) suggests that a state analogous to ~ turbulence ¥ 
exists as a result of creation of mesons during the impact of a particle 
on a nucleus. This would result in an increase of the “ viscosity Fs of 
nuclear matter by which a higher probability should exist for dragging 
out “lumps ” of nucleons from the parent nucleus. a 
In order to determine if there is any correlation between the production 
of energetic heavy particles and meson showers, the stars have been 
divided into different groups according to the values of N,, and number 
of shower particles n,. For each group the ratio between the number of 
stars with energetic helium particles Ny, and the total number of stars 
in the HA. 21 plates N,,4; has been computed. The result is shown in 
Table II. (The absolute values in the Table are in arbitrary units.) 


TABLE II. 


= 


The ratio BES. : 


total 
Np 3-9 10-16 17-23 224 
Ns 
1 1-2+0-8 5-8 +3-6 2144 30+11 
2-3 0-8 -+0-3 4-441-1 2145 50+17 
>4 1-2+0-7 4-641-7 11-43 42412 


There is thus no significant difference in the ratios for groups of stars 
with the same N, and various number of shower particles which indicates 
that no very strong connection exists between the ejection of energetic 
heavy particles and meson showers. The majority of heavy particles 
considered are helium nuclei, but the same conclusions appear to be 


valid for heavier fragments, in accordance with previous results (Perkins 
1950 a). 


(c) Single energetic helium particles. 


A striking feature of the present investigation is the observation of 
single energetic helium particles, starting and ending in the emulsion 
without any associated star. One example has been found with a range 
of 2200 » corresponding to an energy ~ 85 MeV., and four with ranges 
~ 1500 » corresponding to energies ~ 70 MeV. Further single heavy 
tracks are observed to start in the emulsion and pass out of the surfaces. 

Two-prong stars have been observed consisting of an energetic helium 
particle and a short track representing the recoil nucleus of range ~ 5-10 fie 
or with a “thin” track due to a singly charged particle with ionization 
near the minimum value. The latter is presumably the primary, and in 


most cases makes a very small angle with the direction of the heavy 
particle. 


- 
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The single tracks may be either true one-prong stars, or stars in which 
the remaining particles emitted are neutrons. 


Note added in proof.—Recently a two prong star has been found consisting 
of a “ grey” track and the track of a helium particle of range 12,370 p; 
this corresponds to an «-particle energy of 236 MeV. The grain density 
and scattering of the “ grey ” track is consistent with that of a relativistic. 
a-particle, presumably the primary. The angle betweeri the direction of 
motion of the two particles is 16°. 


Discussion. 


Firstly it seems safe to exclude any explanation of the heavy fragments 
based on the nuclear evaporation processes. We are here observing 
particles of kinetic energy far in excess of that due to electrostatic 
repulsion. Further, the angular distributions suggest a process, quite. 
different from that of evaporation, where the particles are emitted almost 
isotropically with respect to the primary. 

At first sight the phenomena are very suggestive of a direct “ knocking 
on ”’ effect of the heavy particle by the primary or by one of the nucleons 
itself recoiling from the primary. Thus the fragment would leave the- 
nucleus before energy had been statistically shared among the nucleons, 
that is in a time very short compared with the “ cooling down ”’ time. 

We select out of the present material all the stars with one or more- 
shower particles. In general primary particles in these stars will have- 
energies E,, in excess of 1000 MeV. Fig. 7 shows a plot of the kinetic 
energy E, of the helium nuclei against y, the angle which they make with 
respect to the primary. The curves A and B represent the function 

4m,,m, 

| =i peeneranie eau Sed eee Se ete Ca) 
calculated from the conservation laws for a free non-relativistic collision 
between the two particles where m, is the mass of an «-particle and m,, 
that of a proton. The curves A and B refer to E,=200 MeV. and 
500 MeV. respectively. It is clear that the vast majority of the x-particles 
cannot have been ejected as a result of a direct collision with the primary,. 
even supposing that both the «-particles and the primary have suffered 
some energy loss in their passage through the nucleus. Thus we conclude: 
that most of the helium fragments have been ejected by encounters with 
secondary nucleons. This conclusion is supported by the observation of 
a-particles at angles greater than 90° to the primary, which must necessarily 


result from secondary effects. 
A similar analysis applied to the lithium fragments gives the same: 


result. ? f 
The energy which an energetic secondary can transfer to the ragmen 
would depend mainly on the relative impact parameter for the collision. 
This is in agreement with the absence of correlation between the energy 


of the heavy particle and the size of the associated star. 
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The close correspondence between the energy spectra: for helium 
particles, protons and deuterons might imply a similarity in their 
‘respective ejection mechanisms. The process by which one nucleon 
can eject an energetic aggregate of nucleons from the parent nucleus 
‘remains to be explained. This phenomenon could perhaps be interpreted 
in terms of long range nuclear forces. | 

A point worth mentioning concerns the masses of the heavy particles 
in the present investigation. It may well be that a large proportion of 
‘the helium particles are of mass 3, since a surprisingly large number of 
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-A plot of the kinetic energy against angle with primary of a-particles from 
stars with n,>1. (For explanation of curves, see text). 


‘energetic tritons are observed to be ejected from stars (Camerini and 
Fowler, private communication). Measurements are now being made to 
‘determine the He/He*-ratio by the scattering method. 

The present observations imply the existence of aggregates of nucleons 
-as temporary sub-units of the nucleus. More detailed investigations may 
possibly throw light on the much discussed problem of the a-particle 
clustering in nuclei (von Weizacker 1938, Wergeland 1941, Wheeler 1941) 


See Ge 
O. C. SORENSEN Phil. Mag. Ser. 7, Vol. 42, Pl. VII. 


A nucleus of Magnesium (Z=12-++1) strikes a silver or bromine nucleus in the 


——F 


emulsion. From the disintegration emerges a slow particle of charge Z=5-+-1. 
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Simultaneous ejection of an energetic helium particle and a lithium fragment 
from a nuclear explosion. 
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XVIII. A Note on Hildebrand’s Approximation for Thermal Pressures 
in Solids. 
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ABSTRACT. 


The validity of a widely used approximation for thermal pressures 
due originally to Hildebrand is discussed. The distinction between the 
“« vibrational ’’ pressure and the thermal pressure is stressed ; Hildebrand ’s 
expression is shown to lie intermediate between the two for all temper- 
atures. Since the percentage difference between the two pressures 
diminishes monotonically with increasing temperature, the accuracy of 
the approximate expression improves with increasing temperature. But 
at temperatures comparable with the Debye temperature, Hildebrand’s 
expression is much closer to the vibrational (10 per cent error) than the 
thermal pressure (40 per cent error). It follows further from the dis- 
cussion that the linearly extrapolated volume from any temperature to 
the absolute zero lies intermediate between the real volume at absolute 
zero of temperature and the volume corresponding to static equilibrium. 


Ir is the purpose of the present note to discuss the validity of an 
approximation, which was first introduced by Hildebrand (1931) and has 
since been widely used in connection with the accurate calculations of 
cohesive energies and interatomic forces of solids (Born and Mayer 1932, 
Mayer and Helmholz 1932, Huggins and Mayer 1933, Fowler 1936, 
Huggins 1937, 1947). It is based on the assumption that the energy of a 
solid can be approximately considered as consisting of a purely volume- 
dependent and a purely temperature-dependent part. Hildebrand’s 
procedure is the following : Let us write the internal energy of a solid as 


EK=U-+E,,,, . . . . . . . ° (1) 
where E,,, and U refer respectively to the vibrational energy and the 


remaining, purely volume-dependent static energy. It is assumed that 
one can approximately write 


dU (dE oF os 
w=(sv),~(sr), *7(Gv), 2 2 
F, 8 being respectively the Helmholtz free energy and the entropy ; 
or, since for normal pressures one can put p= —(0F/AV),—0, 
dU ~T 0s 
Tlie): nD 


T(dS/0V)q is thus used approximately as the “ thermal pressure ”. 


* Communicated by Professor H. Frohlich. 
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Despite the wide usage of the approximation, no pertinent discussions 
of its validity have, however, been offered. Hildebrand (1931) has 
attempted a partial justification empirically, by trying to show that 
alternative determinations of the interatomic forces using the above 
relation at different temperatures, yield the same result; the results 
thus obtained in fact show a distinct drift with the corresponding 
temperatures used. Sometimes it has been implied that (3) is valid 
owing to the small magnitude of the vibrational energy E,,,. This has 
led to the suggestion that (3) is a good approximation at not too high 
temperatures (Fowler 1936) ; the reverse is in fact the case (see below). 

The Mie—Griineisen equation (Griineisen 1926) : 


dU yEvip 4 
Vue site ame Si ian emcee) 


which may, for instance, be deduced by following either the Einstein 
or Debye treatment of vibrations, is ideally suited to a discussion of the 
above approximation. Differentiating (4) with respect to T, we obtain 


(=F ete cin (5) 
1d a 


Moreover, for normal pressures, we can put p=0 in (4): 


dU yE yin 


7 a RY ee 
dN VY HER eye (8) 


OS\ _ (ep ) 
and comparing (3) and (5), we see that the Hildebrand approximation 
for the thermal pressure is given by yT'Cy/V, whereas it follows from (6) that 


yE,;,/V represents the true value; the ratio of the approximate to the 
true value is thus 


Remembering that 


ies 
oa Evin 


Se On wT Ga) 


which, we remark, is equal to unity only for temperatures high enough 
for the classical energy partition law to hold. Nae, Bi ee 

One notes that, owing to the fact that the energy of zeto-pomt vibrations 
depends only on the volume, the “thermal pressure“ can 1n fact be 
defined in two different senses : 7 

(i.) One can include the zero-point energy as part of U. All the above 
formule remain correct as they stand, if E,;, refers only to eat 
vibrations. The thermal pressure defined in this way vanishes at 
T—0° K. (see (6)). 

. P2 
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(ii.) Alternatively, one may, as we have already assumed above, 
take U as the purely static energy. Then the zero-point energy has to 
be included in E,,,. It follows that E,,, approaches a finite value for 
'T=0° K. in this case, and the corresponding ‘‘ thermal pressure ”’ does 
not vanish ; we shall thus refer to it as the vibrational pressure. 

Hildebrand’s approximation is ambiguous on this point, 2. e. whether 
‘T(dS/AV)» should refer to the thermal or vibrational pressure; it is 
probably understood differently by different authors using it. To be 
consistent with the idea of dividing the energy into a volume-dependent 
and temperature-dependent part it would indeed appear the more 
appropriate to identify Hildebrand’s expression T(0S/0V), with the 
thermal pressure rather than the vibrational pressure. For one would 
attempt to include all the purely volume-dependent terms in U, hence 


also the zero-point energy ; moreover, as the approximate expression — 


T(0S/0V), vanishes at T=0°K., it points once more to (i.) as the more 
fitting interpretation. 

Using the Debye theory, one finds that the ratio p defined in (7) is 
only a function of T/O (O the Debye temperature), whether E,,, refers 
to the total vibrational energy or only the thermal part. The curves 
p versus T/@ for both of these alternatives are given in the adjoining 
figure. Curve (i.) corresponds to the case that E,,, refers only to the 
‘thermal vibrations, and curve (ii.) refers to the alternative case, where 
Ey, is the total vibrational energy. If we interpret Hildebrand’s 
‘expression as referring to the thermal pressure (see (i.)), curve (i.) gives 
the ratio of the approximate to the true value ; if, on the other hand, 
we interpret Hildebrand’s expression as referring to the vibrational 
‘pressure (see (ii.)), then curve (ii.) represents the proper ratio. Both curves 
are seen to approach the line of unity towards high temperatures ; 
-however, for the temperature range of practical interest, the Hildebrand’s 
‘expression is much more accurate, if interpreted as the vibrational 
‘pressure (curve (ii.)) rather than the purely thermal pressure (curve (i.)). 
This is exactly the reverse of the more naive anticipation explained earlier. 
For all the alkali-halides except LiF and NaF, Hildebrand’s expression 
used at room-temperatures will be within 10 per cent of the actual 
vibrational pressure. 

An interesting conclusion can be drawn from considering the equilibrium 


V,) defined by the equation 
7) 
(a Ms =—() . . . . . . . . (8) 


4 
with the approximate equation (3). By expanding dU/dV in (3) with 
respect to V—V, about the above equilibrium position and neglecting 
second and rates order terms, we find peti: that 


T(S/8V) 


Vo2V— (ZU/dv?), . . . . . . . . (9) 


— 


a 


’ 
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The interpretation of Vy from the point of view of this approximation 
is as ambiguous as the interpretation of Hildebrand’s expression itself : 
if the latter is taken as the thermal pressure (i.), then U should include 
the zero point energy and the corresponding V, is the volume of the 
solid at T=0° K; but if it is taken as the vibrational pressure (ii.), U is 
the static energy and V, is the static equilibrium volume corresponding 
to the hypothetical case that the atoms were at rest. The fact that the 
curves (i.) and (ii.) lie respectively above and below the line of unity shows 
that Hildebrand’s expression has a value intermediate between the true- 
thermal pressure and the true vibrational pressure, whatever the 
temperature. It follows, as a moment’s consideration will show, that. 


e 


0-5 T/@ | IBS} 2 
Ratio of Hildebrand’s pressure to the thermal (i.) and. vibrational pressure (ii.). 


the value of V, as given by (9) lies in fact between the equilibrium 
volume of the static case and the actual volume at the absolute zero of 
temperature. Moreover, if the relation (9) is used at room seh aed 
V, should be much nearer to the static equilibrium value than the va a 
at T—0°K, for we have seen that the approximate expression shou 
be much nearer to the vibrational than the thermal pressure. 

Now if we write on the right-hand side of (9) 


(=a), ~ (Bh), = “ort = ~ (>0), (5), 
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and put approximately 


Op e? : 
(3), lan (Fr), , (difference ~10 per cent) 


(9) becomes simply 


av 
VeeV—1 (55) . oo) Lae 


This means that the equilibrium volume V, deduced from Hildebrand’s 
approximation is practically identical with that obtained by a simple 
linear extrapolation to the absolute zero of temperature. Hence we can 
express the conclusion of the last paragraph in the rather interesting 
form: the linearly extrapolated value for the volume of a solid from any 
temperature to the absolute zero lies in between the static equilibrium volume 
and the actual volume at absolute zero of temperature. Moreover, if it is 
extrapolated from a temperature comparable or higher than the Debye 
temperature of the solid, the result is expected to lie much nearer to the 
static equilibrium value. 


REFERENCES. 


Born, M., and Mayzr, J. E., 1932, Zeit f. Phys., 75, 1. 

Fowter, R. H., 1936, Statistical Mechanics (Cambridge: University Press), 
2nd Ed., p. 324. 

GRUNEISEN, E., 1926, Handb. d. Phys., 10, 1. 

HILDEBRAND, J. H., 1931, Zeit. f. Phys., 67, 127. 

Hueerns, M. L., and Mayr, J. E., 1933, J. Chem. Phys., 1, 643. 

Houcerns, M. L., 1937, J. Chem. Phys., 5, 143 ; 1947, J. Chem. Phys., 15, 212. 

Mayer, J. E., and Hetmuotz, L., 1932, Zeit. f. Phys., 75, 19. 


[ 207: ] 
XIX. The Emission of Li8 from Cosmic Ray induced Nuclear Disintegrations. 
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SUMMARY. 


Lithium 8 nuclei ejected from cosmic-ray induced disintegrations have 
been investigated using the photographic plate method. Earlier work 
has been confirmed, and the frequency and energy distribution of the Li® 
nuclei determined and compared with the predictions of the evaporation 
theory. Some nuclear reactions involving the production of Li® have 
been identified. 


§ 1. INTRODUCTION. 


» In 1947, Occhialini and Powell observed a T-shaped track caused by 
a particle emitted from a nuclear disintegration in a photographic plate 
exposed to cosmic rays. These events are also known as “ hammer ” 
tracks. 

Franzinetti and Payne (1948) analysed thirty examples and confirmed 
Occhialini and Powell’s conclusion that they are due to the emission of 
lithium 8 nuclei which decay into two alpha-particles and an electron, 
although the electron track was not observed in their less sensitive plates. 
This was later observed by Occhialini and by Pickup and Voyvodie (1949) 
in emulsions sensitive to particles of minimum ionization. 

This work has been confirmed, and, in addition, the frequency and 
energy distributions of the Li® particles determined. 


§ 2, EXPERIMENTAL PROCEDURE AND SUMMARY OF RESULTS. 


Several batches of Ilford Nuclear research plates, coated with C2 and 
G5 emulsions, and exposed to the cosmic rays on the Jungfraujoch, have 
been examined for hammer tracks associated with stars. 

The energy distribution of the secondary alpha-particles was found 
from their ranges in the emulsion and it agreed with the results of 

'Franzinetti and Payne (1948) for hammer tracks associated with cosmic 
tay stars and those of Baxter, Burcham and Paul (1950) and Rumbaugh, 
Roberts and Hafstad (1937) for Li® nuclei produced artificially. 

Electron tracks were found associated with most of the hammer tracks 
in electron-sensitive emulsions. 

The energy distribution of the ejected particles was measured and is 

i in § 4. 
ae distribution of the hammer tracks was found to be 
consistent with isotropic emission of the Li® nuclei. 


eel 


* Communicated by Sir George Thomson. 
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The frequency of emission of particles forming hammer tracks is given 
in §5 and compared with the predictions of the evaporation theory for 
lithium 8 nuclei of zero spin. 


§ 3. NATURE OF THE PARTICLE CAUSING THE HAMMER TRACKS. 


Tracks of decay electrons were found associated with most of the 
hammer tracks occurring in electron-sensitive emulsions. The electron is 
energetic enough to be at or near minimum ionization and yet slow enough 
to be quite considerably scattered. The electron tracks are therefore 
difficult to see and are in general only found if the emulsion is developed 
to its fullest sensitivity. 

Hammer tracks associated with a decay electron may be due either to 
Lif or, as recently suggested by Alvarez (1950), B5 :— 


Li > Be§*+e°,, Be&* + 2Hef, 
BS > Be§f*+e°.,, Be§* — 2Heé. 


Because of the higher charge of B3 its emission probability will be less 
than that of LiS. 

The experimental results are not incompatible with some of the hammer 
tracks being without a decay electron, and attributable to Bef. In this 
case the excited state of Be? must be a metastable one of long life since 
the life-time of the broad state at 3 MeV. is far too short to allow the 
nucleus to come to rest in the emulsion before disintegrating. 

Since Li§ and B3 decay via the same excited state of Be§ (Alvarez 1950) 
no distinction can be made between them from the energy distribution of 
the secondary alpha-particles. 

It seems likely that the majority of the particles producing hammer 
tracks are Li§ and in what follows it will be assumed that they are all 
attributable to this nucleus. 


§4. THe Enercy DistriputTion oF THE Li8 NUCLEI. 


The energies of the lithium nuclei were determined from their ranges 
and a range-energy curve derived from that of Bradner, Smith, Barkas 
and Bishop (1950) for protons by the relation 


vont) 


where KE, R, « and m are the energy, residual range, charge and mass of 
the particle. At very low energies the lithium nucleus picks up electrons 
and reduces its effective charge. The above formula cannot be used in 
this region with accuracy. But, since this lowering of effective charge 
takes place, for lithium nuclei, only over the last two or three microns, 
its effect is small and has been neglected. 

It was found that the average energy of the lithium nuclei from large 
stars was appreciably greater than that of those from small stars. 
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This may be explained by assuming that the lithium nuclei from large 
stars are ejected from energetic disintegrations of silver or bromine, while 
the lithium nuclei from small stars are the residues after the ejection of 
protons and alpha-particles in the disintegration of carbon, oxygen or 
nitrogen. This interpretation has also been advanced recently by 
Courtenay Wright (1950) after examination of Li§ nuclei associated 
with stars produced in various gases by the cyclotron beam. 

In order to test this hypothesis, all stars associated with hammer tracks 
were separated into two classes, namely those with more than four tracks 
(including the hammer track), and those with less. The additional 


Fig. 1. 
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Energy of Li® in MeV: 
Energy Distribution of Lithium nuclei. 


condition was applied that stars with a recoil, a short heavily ionizing 
track characteristic of the disintegration of a heavy nucleus, were included 
in the former class, irrespective of their numbers of tracks. These classes 
are called A and B. In addition, it is possible for a lithium pucioue to 
be the residue after the complete disintegration of a ia or pos 
nucleus. Only one example has been found, in a star of 27 tracks, and it 
i in class B. i ¢ 
eee se sibistiok of the lithium nuclei is shown in sas sé ane two 
classes being distinguished by shading. The observed distribution was 
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corrected for those long tracks which were lost because they passed out 
of the emulsion. The disparity between the energy distributions of 
lithium nuclei in the two classes is strong support for the above hypothesis 
on the mechanism of their production. 


§ 5. FREQUENCY OF EMISSION. 


In all, it was found that 74 hammer tracks were associated with stars 
in emulsion containing 23,700 stars of more than two tracks. Hammer 
tracks associated with stars caused by 7~-mesons are omitted from these 
figures and will be considered later (§ 7). The total emission frequencies, 
for C2 and G5 emulsions, are showm in Table I. The value for G5 
emulsion is less than that for C2 on account of the greater number of 
stars with more than two tracks. 

The fraction of hammer tracks in class A is also given in Table I. 


TABLE I. 


Observed Corrected 


Hammer tracks—No. of stars >2 tracks | 0-0032-L0-0007 | 0-0052--0-0012 
in C2 emulsion 


Hammer tracks— No. of stars >2 tracks | 0-0023-+-0-0005 | 0-0034-L0-0006 
in G5 emulsion 
Fraction of hammer tracks of class A 0:27 +0-07 0-33 +0-08 


The frequency of emission of particles causing hammer tracks is plotted 
in fig. 2 as a function of star size and compared with the calculations on 
the evaporation theory by Le Couteur for the emission of lithium 8 of 
zero spin from nuclear disintegrations of various sizes. Comparison of 
the observed and calculated frequencies, which should only be made 
for hammer tracks of class A, shows good agreement. 

Hammer tracks are lost if they pass out of the emulsion (see § 4) and 
also if the angle between the projected directions of the lithium and the 
two alpha-particles is too small to be resolved. The frequencies of 
emission have been corrected for tracks lost in both these ways. 

From the frequencies of emission given above, and assuming an 
approximate value of the intensity of star-producing radiation at the 
Jungfraujoch, values of the cross-section for production of Li§ from light 
nuclei can be deduced which are of the same order of magnitude as those 
given by Courtenay Wright (1950) and Titterton (1950). 


§ 6. NUCLEAR REACTIONS IN THE LIGHT ELEMENTS. 


On the assumption that the hammer tracks of class B come from the 
disintegration of light nuclei it is often possible, when the event occurs 
in C2 emulsion, to identify the nuclear reaction concerned. The identifi- 
cation is subject to some uncertainty because lightly ionizing particles 


en ee 
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are not visible in C2 emulsion, and protons, deuterons and tritons are 
not always separable. Analysis of events in G5 emulsion was not 
attempted because of the difficulty of distinguishing between protons and 
alpha-particles. It is possible to confirm the reactions by measurement 
of energy and momentum balance, providing no neutrons are emitted. 
Fig. 2. 
0-1 
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Number of evaporated particles. 
Fraction of stars with hammer tracks as a function of star size. 


The following reactions have been identified 


n+C? > Li§+He+H} (Six examples) 
nu +C? > Li§+3H}+2n9 
or ni-+Cl2 > Li8+2H}+H}-+n) (Eight examples) 


or ni-+Cl2 > Li§+-H}+ 2H} 
ni+Ni! > Li8+ He+2Hj +75 

or ni-+N}4 — Li8+He$+Hj+H; 
ni+Ol > Lif+2He5+ Hj (One example) 


| (Two examples) 
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§ 7. Liratum NUCLEI EMITTED FROM STARS INITIATED BY a—-MESONS. 


In addition to the hammer tracks analysed above, six have been 
observed associated with stars produced by 7 ~-mesons. The total 
number of such mesons in the emulsion examined was 1160. The fraction 
of meson-produced stars with hammer tracks is therefor 0-005 +0-008, 
a value which agrees with the 0-003 given by Bradner (1950) and the 
11/3000 or 0:004-+40-001 found by Adelman and Jones (1949). 

The following nuclear reactions have been identified :— 


Ol +79 , > Li§+He§+2H}+ 2nj 
or O16-+-7° , > Li§-+He§+ H!+H?+-nj 
CP+7r®, > Li§+2H}+2n} 
or C247), — Li§+Hi+H7+nj 


3 (Two examples) 


(One example) 
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XX. The Law of Constant Resolved Shear Stress in Crystal Plasticity. 


By F. R. N. Nasarro. 
Department of Metallurgy, The University of Birmingham*, 


[Received December 27, 1950.] 


ABSTRACT. 


The formula of Peach and Koehler for the total force acting on a 
dislocation in a stressed crystal is derived by the use of a theorem due to 
Colonnetti. The same method distinguishes conservative motions of the 
dislocation, which do not alter the volume, from other motions. The 
component of the total force which produces conservative motions acts 
in the glide plane, and is proportional to the same component of the total 
Stress as appears in the law of resolved shear stress discovered 
experimentally by Schmid. 


Tue forces acting on dislocations are most conveniently derived by the 
use of a theorem due to Colonnetti (1915). Consider a body containing 
a closed loop of dislocation with Burgers vector b, and cover this loop 
with a cap. Let an element of this cap have area dz, and let the body 
be subject to external stresses as a result of which the force acting across 
this element d= of the cap is Xd. The signs of the vectors b, d& and 
_ do must be determined conventionally. The convention used here is 
that of Peach and Koehler (1950). Colonnetti’s theorem then states 
that the work W, done by the external forces when the dislocation is 
formed is given by 


Vy ea cage caeeee os te (1) 
If the stress produced by the external forces is 38, X is given by 
Rea eee ro. sr ee we (2) 
and 
Weenie eee tee te (9) 


Suppose now that an element do of the dislocation moves a distance dx. 
The area of the dislocation loop is increased by 
62=dx a do, mated Bderriy a0) da, ig (4) 
and the external forces do work 
dW ,=b . (dx ado) 


=[(by)ado].dx.. . . . . + + (5) 
‘The force dF on the element de is thus given by 
GED do eee tit veel: 4 (8) 


This is the result of Peach and Koehler. 
* Communicated by the Author. 
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It is now necessary to distinguish between those motions of the 
dislocation which conserve the volume of the crystal (conservative 
motions) and those motions which produce interstitial ions or vacancies 
(non-conservative motions). Only conservative motions can lead to 
steady slip, though non-conservative motions may be important in 
work hardening and especially in thermal recovery. The necessary and 
sufficient condition that a motion should be conservative is that a uniform 
hydrostatic pressure p should do no work when the motion occurs. For 
such a pressure, b39 is given by b33=bp, and it follows from (5) that the 
motion is conservative if 


(ba do)iidx=20.) 0 eo atl Seay) 


Thus the motion is always conservative if b~ado=0, that is to say, if 
the element of dislocation is of pure screw type. If b and deo are not 
parallel, they define a plane, the glide plane, and the motion is conservative 
only if dx lies in this plane. The normal to the glide plane is the unit. 
vector n defined by 


n=(bado)/|ba dol. e. 0): -« a nee eee (8) 


The component df of the force dF which is effective in producing slip 
is the component lying in the glide plane, given by 


df=d¥F —n(n dF) 
ssxet a (ma @B).9 )-.imorss 5.08 eee ee 
Now * 
na dF=n a [(b) a do] 


=(bs)(n do)—do[n(bP)] 


=—doa(bPn), 
and so 


df=(n,do)(bBn). . . ... «. « « « (10) 


Since $n is the force acting across unit area of the glide plane, the force 
per unit length tending to produce glide is 6 times the component in the 
glide direction of the force acting across unit area of the glide plane, and 
acts in the glide plane normal to the line of the dislocation. If a loop 
of dislocation lies entirely in the glide plane, its conservative motions are 
confined to this plane (except for those elements which are of pure screw 
type), and its whole length is subjected to the same normal force tending 
to increase the area of the loop. This is in accordance with the experi- 
mental result (Schmid 1924) that slip begins when the “ resolved shear 
stress ’’ bj3n/b reaches a value which is characteristic of the crystal, but 
independent of the other components of the applied stress. 
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XXI. CORRESPONDENCE 


The Production of Neutrons by High Energy Protons. 


By J. M. Casszts, T. C. Ranpuz, T. G. Pickavance and A. E. TaYLor*, 
Atomic Energy Research Establishment, Harwell. 


[Received December 8, 1950.] 


THE production of neutrons by high energy protons has been investigated 
with the help of a differential detector. The third counter in the triple 
coincidence telescope used previously for total cross-section work (Taylor, 
Pickavance, Cassels and Randle 1950) was replaced by a gridded ionization 
chamber sensitive only to comparatively slow protons. It followed that. 
protons ejected from the polythene radiator were counted only if they 
possessed energies lying between certain limits set by a carbon absorber 
placed immediately in front of the chamber. The counts resulting from 
n, p collisions in the hydrogen content of the polythene were separated out 
by background runs with a carbon radiator. After correcting for protons 
lost while slowing down in the absorber, absolute neutron fluxes in the 
range 50-200 MeV. could be found with the help of the Berkeley 
experiments on the n,p interaction (Hadley, Kelly, Leith, Segré, 
Wiegand and York 1949, Brueckner, Hartsough, Hayward and Powell 
1949, Kelly, Leith, Segré, and Wiegand 1950). 

The correct operation of the telescope was checked by using it to deter- 
mine the total cross sections of carbon and hydrogen at 156 and 110 MeV. 
The results were consistent with those of other work at Berkeley and at 
Harwell, showing that neutrons of the intended energy were being detected. 

The energy spectra of the neutrons produced by bombarding beryllium, 
carbon, aluminium and uranium targets were then investigated (figs. 1-4). 
These targets had thicknesses of 1-67, 1-09, 1:19 and 1-82 grms.-cm.~* 
respectively. The resolution diagrams in fig. 1 show the ‘“‘smearing ” effect 
of the neutron detector only. There is additional lack of resolution 
caused by the slowing down in the target of the primary protons, whose 
energy was initially 171 MeV. The magnitude of this can be gauged from 
the mean proton energies set out in the Table, column 2. No corrections for 
‘‘smearing’’ have been applied to the results given. The fluxes shown in 
the diagrams refer to a point 15-3 metres from the target at an angle of 23° 
from the forward direction. They are normalized to an actual current of 
1 microampere being accelerated by the cyclotron, but the effective 
current was greater than this because protons could traverse the target 
several times before being lost. It was found to be 3-4 microamperes in the 
carbon target by measurement of the rate of production of 1C, a process 


* Communicated by the Authors. 
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Figs. 1-2. Neutron spectra obtained from beryllium, carbon, aluminium 
and uranium by bombardment with 171 MeV. protons. 
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' Fig. 3-4. Neutron spectra obtained from beryllium, carbon, aluminium 
and uranium by bombardment with 171 MeV. protons. 
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for which the cross-section is known from work at Berkeley (Aamodt, 
Peterson and Phillips, private communication, 27th February, 1950). 
Knowing this effective current, it was possible to calculate the differential 
cross-section, at anangle of 24° from the forward direction, for the production 
in carbon of neutrons with energy greater than 60 MeV. (see Table, column 3). 
If a constant fraction of the cyclotron beam was lost at each traversal 
of the target, then the mean energy of the primary protons was 155 MeV. 

The loss of protons in the target is determined by multiple Coulomb 
scattering (Moyer, private communication), and it has been found for the 
Harwell cyclotron that the number of times the proton beam passes through 
a target at a given orbit radius is such as to bring the R.M.S. angle of 
multiple scattering up to a constant value, provided that the target is 
thin enough to make several traversals necessary. Using this principle, 
the neutron production cross-sections in beryllium and aluminium were 
obtained. The mean proton energies in these two targets were 156 


Mean proton energy | Production cross-sections x 107 . 


Target element (MeV.) (cm?) 
24° ate 
Beryllium 156 91+ 45 99+19 
— 30 
Carbon 155 47+. 23 72+15 
— 16 
Aluminium 166 122+ 61 — 
— 41 
Uranium 166 1520+760. _ 
—507 


and 166 MeV respectively. The uranium target was thick enough to 
suggest that the beam went through it only once, and the cross-section and 
mean proton energy given were calculated on this basis. These four cross- 
sections may be in error by a factor of 1-5 either way because systematic 
errors could have been introduced in estimating the cyclotron beam 
current during the various runs. 

The neutron production cross-sections in beryllium and carbon targets 
at an angle of 5° to the forward direction were also measured, in another 
and more direct way (see Table, column 4). The gridded chamber in the 
telescope was replaced by a counter sensitive to high energy protons and 
the carbon absorber adjusted so that all neutrons over 60 MeV. were 
counted. The beryllium target for this experiment consisted of two thick 
sheets of beryllium interleaved with three thin sheets of carbon. After 
two short runs in which neutrons produced were counted for each target in 
turn, the targets were removed from the cyclotron and the ™C content of 
each measured with a calibrated counter. The background measurements 
required to find the fraction of counts in the telescope due to hydrogen 
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recoils had been done on the previous day, so that effectively no C counts 
remained from those runs at the time of the calibration run. With the 
help of the known ”C(p, pn)!!C cross-section, the neutron production 
cross-sections could be found. The stated standard error is intended to 
include all systematic errors. 

The shape of the beryllium spectrum reported here seems quite consistent 
with the results of other experiments using protons of 110 MeV. (Bodansky 
1950) and 340 MeV. (Fox, Leith, Wouters and MacKenzie 1950). 

We hope to improve and extend these results and publish them, with a 
more detailed account, at some future time. 
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The Velocity Brownian Movement. 


By D. K. C. MacDonatp *, 
Clarendon Laboratory, Oxford. 


[Received December 14, 1950.] 


Dr. Firru’s amplifying and detailed discussion (Fiirth 1950) of the 
example I appended to illustrate the main argument of my paper 
(MacDonald 1950) is welcome, but does not, of course, alter any of the 
fundamental conclusions. Nevertheless, may | make the following brief 
comments :— 

(1) ‘The purpose of my example, which it appears Dr. Firth may have 
misunderstood, was to show that the available statistical data on the 
observed velocity of a Brownian particle would indeed be already sufficient 
to determine D (and hence k=677aD/T), at least in order of magnitude, if a 
reasonable value for the discrimination-period of human vision were 
assumed. Since I was in fact only concerned with demonstrating the 
underlying principle I deliberately chose for simplicity the familiar one- 
dimensional Brownian movement equation. 

(2) The use of the standard deviation to apPron ac the mean 
absolute velocity is clearly a good approximation for a’ stationary system 
in thermal equilibrium as for the Brownian particle (or in the general 


* Communicated by the Author, . 
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kinetic theory of gases) but, of course, for a system with an externally 
superposed mean velocity (as in a current-carrying medium) the two need 
have no connection with one another. 

(3) The thesis underlying Dr. Fiirth’s first two sections is not entirely 
clear to me. However, no physical variable can have a frequency- 
spectrum extending without limit and consequently, of course, can never be 
perfectly random (and therefore non-differentiable). It is nevertheless 
quite customary to refer to a quantity as varying stochastically if within 
the frequency-range of observation the “ power ”’ spectrum is effectively 
flat as, for example, in the case of thermal noise from a resistor at all 
accessible radio frequencies and normal temperatures. 
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Structural Chemistry of Inorganic Compounds, Vol. I. By W. Htoxet. 
(Elsevier Publishing Co. Inc.) Price $9. (£3 10s. Od.) 


Tuts is essentially a book for use in the latter stages of a University degree 
course. About one-eighth of the book is taken up by an excellent historical 
survey. The remainder is about equally divided between a presentation of the 
coordination theory which is the basis of the author’s systematization, and a 
discussion of atomic structure and physical methods for investigating molecular 
structure. This order is a little unusual, and emphasizes the importance 
attached to stereo-chemistry and the coordination approach. 

The author can be congratulated on the large amount of recent information 
which he has included, though no mention is made of the optical resolution of the 
dodecammine-hexol-tetracobaltic salts, the first resolution of a purely inorganic 
compound. An incorrect structure is assigned to this group of salts. The 
correct structure has in fact been determined by X-ray analysis, contrary to the 
statement on p. 166. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. 
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